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The  SCF-Xa  method  has  been  applied  to  the  study  of  crystalline  defects. 
First  the  LiF  molecule  was  studied.  The  Xa  method  gave  the  correct  inter- 
nuclear  separation  and  vibrational  frequency.  This  work  was  extended  to 
the  LiF  crystal  via  a cluster  treatment.  The  correct  lattice  constant 
was  derived  and  the  electronic  structure  was  in  agreement  with  experiment. 
The  final  phase  of  the  work  dealt  with  the  removal  of  an  F"  ion  creating 
an  F - center.  The  Xa  - cluster  method  gave  a good  description  of  this 


FORM 

I JAN  n 


EDITION  OF  I NOV  EE  IS  OBSOLETE 
S/N  0102-014- E60I 


SECURITY  C 


ibXO  0 


ilFICATION  OF  THIS  RAOE  fWh»n 


iimup.  i‘  ■•■■■IT- — ^ — — ^ — » 


- 1 - 


ONR  FINAL  REPORT 


This  is  the  Final  Report  of  ONR  Contract  N00014-75-C-0785.  The  main 
focus  of  the  research  has  been  on  solid  state  theory  with  an  emphasis  on  the 
magnetic  properties  of  matter  and  defect  in  crystals.  Magnetic  phase  transi- 
tions, energy  hand  structures,  and  the  formation  of  vacancies  and  interstitials 
were  studied.  In  addition  very  large  magnetic  pulses  suitable  for  remote 
electromagnetic  detection  were  generated  in  a joint  program  with  the  National 
Magnet  Laboratory.  These  pulses  had  peak  currents  of  22,000  ampues,  pulse 
lengths  from  0.1  uc  to  30  usee,  and  maximum  pulse  power  of  48  mega  watts.  The 
following  Ph.D.  theses  were  completed  under  this  contract. 

1.  Eric  Li,  "Electronic  Energy  Structure  of  FgS2  and  N^S2" 

2.  D.  Choo,  "SCF  Xa  Cluster  Model  Calculations  for  Crystals 

and  Point  Defects" 

A detailed  study  of  the  electronic  states  associated  with  impurities, 
iso-electronic  traps,  and  vacancy  and  interstial  defects  was  made.  A cluster 
technique  developed  by  Professor  Keith  Johnson  was  extended  and  applied  to 
these  problems.  This  work  is  included  in  this  report  as  Appendices  I - III. 

The  objective  of  this  program  was  to  develop  a technique  adequate  in  scope 
to  handle  defect  problems  in  nuclear  reactor  materials  where  the  condensation 
of  individual  lattice  vacancies  into  voids  and  then  a lattice  of  voids  leads 
to  a swelling  of  materials  exposed  to  high  levels  of  radiation.  This  swelling, 
by  as  much  as  25%,  leads  to  a serious  structural  weakening  and  hence  to  a dang- 
erous made  of  operation.  The  ability  of  the  SCF  - Xa  method  to  describe  va- 
cancies and  interstitials  suggests  that  it  may  be  useful  in  studying  radiation 
damage  in  nuclear  reactor  materials. 
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Perhaps,  one  of  the  most  basic  defects  is  the  F - center  in  alholi  bolides. 

LiF  was  selected  as  a practical  case  to  initiate  the  research.  First  the  LiF 

i 

molecule  was  treated,  then  the  LiF  crystal  simulated  by  a Li4f^4  cluster  and  ] 

i 

finally  one  F was  removed  to  study  the  F center.  The  results  obtained  are 
stated  in  detail  in  the  appendices. 

The  principal  investigator  would  like  to  express  his  appreciation  for 
the  support  of  the  ONR  in  this  research.  A proposal  is  now  being  considered 
by  ERDA  to  continue  with  the  study  of  radiation  damage. 
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SCF  Xa  CLUSTER  MODEL  CALCULATIONS  FOR  LiF 
MOLECULE,  LiF  CRYSTAL,  AND  ITS  F-CENTER* 


D.  Choo^  and  G.W.  Pratt,  Jr. 

Department  of  Electrical  Engineering  and  Computer  Science,  and 
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ABSTRACT 

The  SCF  Xa  multiple-scattering  method  of  Johnson 
has  been  applied  to  small,  neutral  cluster  models  of 
a LiF  molecule  and  a LiF  crystal  with  an  overlapping- 
sphere  scheme.  In  both  cases  the  calculated  bulk  pro- 
perties at  equilibrium, such  as  internuclear  distance, 
vibrational  quantum  energy  or  bulk  modulus,  and  cohes- 
ive energy,  are  in  very  good  agreement  wth  experimental 
results.  The  F-center  in  LiF  crystal  is  also  studied 
using  the  same  cluster  model  of  the  crystal. 
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I.  INTRODUCTION 

The  purpose  of  the  present  work^  is  to  model  a crystal- 
line solid  by  a small  cluster  of  atoms  in  order  to  study  the 
electronic  structures  of  the  perfect  crystal  and  the  crystal 

with  a point  defect  via  the  SCF  Xct  multiple-scattering  method 
2-7 

of  Johnson  from  first  principles. 

Over  the  decades  various  energy-band  methods,  such  as 
APW,  OPW,  pseudo-potential,  and  KKR  methods,  have  been  devel- 
oped and  applied  to  many  perfect  crystals  for  their  electronic 
structures  which  includes  long-range  periodic  effects.  Still 
some  basic  equilibrium  properties  of  a crystal  can  be  obtained 
from  a small  portion  of  the  crystal  through  its  short-range 
effects,  which  will  i^e  shown  in  this  paper.  Such  a model  will 
also  be  ideal  for  investigating  point  defects  in  the  crystal. 

It  has  been  possible  to  solve  the  point-defect  problems 

from  first  principles  by  the  energy-band  method  through  the 

0 

Green’s  function  formalism  of  Koster  and  Slater  ,e.g.  the  single 

9 

and  di-vacancies  in  silicon  by  Callaway  and  Huges  , and  the 
Pb  and  Te  vacancies  in  lead  telluride  by  Parada  and  Pratt^*^. 
Calculation  was  also  carried  out  for  the  F-center  in  LiF  by 
the  LCAO  method  recently  by  Chaney  and  Lin^^.  The  SCF  Xci  multiple- 
scattering method  has  great  advantages  over  such  traditional 
methods  in  its  straight-forwardness  and  ease  in  self-consiscent 
computation  with  rapid  convergence.  However,  in  the  applications 
of  the  method  to  crystals,  careful  considerations  must  be  given 


3 


to  the  boundary  conditions  of  the  cluster  models  because  of  the 

inherent  molecular  nature  of  the  clusters.  This  method  has 

been  applied  to  calculate  the  electronic  properties  of  vacancies 

in  PbS,  PbTe,  and  SnTe  using  relatively  large  clusters  of  27 

atoms  in  order  to  make  the.m  nearly  neutral  including  several 
12 

atomic  layers. 

LiF  has  both  molecular  and  crystalline  forms  so  that  the 
LiF  molecule  can  be  first  tested  on  how  much  and  how  well  we 
can  obtain  its  properties  by  the  SCF  Xct  multiple-scattering 
method  before  we  go  into  simulating  the  LiF  crystal  by  a cluster 
model  to  calculate  the  crystal  bulk  properties.  The  changes 
of  the  properties  as  the  LiF-  molecules  are  transformed  into 
crystal  can  be  also  examined.  The  constituent  atoms  of  LiF 
have  relatively  small  n-jmbers  of  electrons,  therefore,  the  cal- 
culations can  be  performed  easily.  Also  alkali  halides  in 
general  have  been  studied  rather  thoroughly,  both  theoretically 
and  experimentally,  so  that  there  are  much  data  already  avail- 
able includi.ng  the  knowledge  about  their  point  defects.  The 
F-center  problem  in  the  LiF  crystal  is  also  studied  using  the 
cluster  model  for  the  perfect  crystal.  The  LiF  is  almost  en- 
tirely ionic  in  its  bonding  character;  more  covalent  crystals 
are  treated  in  the  following  papers.^ 

The  LiF  molecule  is  first  treated  in  Section  II,  and  the 
perfect  LiF  crystal  in  Section  III,  followed  by  the  F-center 
problem  in  the  LiF  crystal  in  Section  IV. 
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II.  LiF  MOLECULE 

A.  Calculational  Considerations 

The  heteronuclear  diatonic  molecule  LiF  is  mostly  ionic 

in  its  bonding  character  when  it  is  in  equilibrium  with  the 

o 13 

internuclear  distance  of  2.354  a.u.  (1.51  A).  In  order  to 
determine  the  geometry  of  the  LiF  molecular  cluster,  we  first 
consider  the  sizes  of  the  Li  and  F atomic  spheres.  The  atomic 
(neutral)  and  ionic  radii  of  Li  and  F which  were  deduced  from 
crystal  data  are  shown  in  Table  I.  The  atomic  radii  are  so 
much  different  from  the  ionic  radii  that  they  will  not  be 
suitable  in  describing  the  LiF  molecule  (or  LiF  crystal  for 
that  matter) . 

Suppose  we  make  the  Li  and  F atomic  spheres  touch  each 
other  to  form  the  molecular  cluster,  keeping  the  raio  of  their 
ionic  radii  constant.  Then  the  radius  of  the  F atomic  sphere 
becomes  1.888  a.u.,  while  the  radius  of  the  Li  atomic  sphere 
becomes  0.965  a.u.,  which  is  too  small  even  for  the  Li"*”  ion. 
Furthermore,  the  interatomic  region  becomes  relatively  too 
large  occupying  about  two  thirds  of  the  total  volume  enclosed 
by  the  outersphere.  In  muffin-tin  approximation  of  potentials, 
which  is  adopted  in  the  SCF  Xa  multiple-scattering  method,  the 
interatomic  potential  is  only  volume-averaged  and  consequently 
less  accurate  than  the  atomic  potentials  which  are  spherically 
averaged.  In  order  to  correct  some  errors  which  arises  from 
this  geometry,  we  employ  an  overlapping  scheme  in  which  the 


spherical  atomic  regions  are  allowed  to  extend  beyond  the 

situation  of  the  touching  spheres,  thereby  overlapping  with 

neighboring  atomic  regions.  Similar  overlapping-sphere 

schemes  have  been  applied  to  molecules  resulting  in  better 

agreement  with  experiments  than  the  cases  without  using  such 
14  15 

schemes.  ' It  will  also  allow  the  atomic  spherical  regions 
to  be  large  enough  to  accommodate  most  of  the  electrons  which 
are  primarily  associated  with  the  atoms.  The  overlapping  is, 
however,  not  to  be  carried  too  far  because,  if  an  atomic 
center  gets  deep  into  the  core  region  of  a neighboring  atomic 
center,  its  core  charge  will  adversely  affect  the  spherically 
averaged  potential  in  the  outer  region  of  the  atomic  center 
under  consideration. 

For  the  LiF  molecular  cluster,  we  start  from  the  ionic 
radii  which  results  in  an  overlapping  scale  factor  1.33  and 
varies  it  such  that  the  virial  theorem  is  satisfied.  It  was 
shown  by  Slater^®  that  the  virial  theorem  is  exactly  satisfied 
by  the  SCF  Xa  multiple-scattering  method.  The  outer-sphere 
is  always  adjusted  to  touch  both  atomic  spheres.  The  over- 
lapping scale  factor  for  the  LiF  molecular  cluster  is  fixed 
at  1.4  after  the  first  few  iterations  at  which  point  the 
virial  theorem  is  satisfied  to  within  0.2%.  If  we  satisfy 
the  virail  theorem  more  accurately,  the  overlapping  scale 
factor  will  be  larger  and  the  core  charge  of  a neighboring  atom 


start  to  affect  the  potential  erroneously  by  a significant 
amount.  The  same  overlapping  scale  factor  is  used  consistently 
for  the  molecular  cluster  with  various  internuclear  distances. 
Fig.  1 shows  the  cluster  geometry  for  the  LiF  molecule. 

The  atomic  exchange  parameter  a of  Schwarz^"^  are  used  for 
the  corresponding  atomic  centers  of  the  cluster  (0.78147  for 
Li,  and  0.73732  for  F) . For  the  interatomic  and  extramolecular 
regions,  a weighted  average  of  the  atomic  a's  with  respect  to 
the  number  of  valence  electrons  in  the  atoms  is  adopted  since 
we  are  primarily  interested  in  the  electron  states  of  the 
valence  electrons.  It  is  0.74468  for  the  LiF  molecule  (and 
the  same  for  the  LiF  cyrstal). 

The  molecular  basis  functions  are  made  of  Li  s,p  and  F 
s,p  atomic  functions.  The  contributions  from  the  excited 
atomic  levels  are  negligible  due  to  the  rapid  convergence  of 
the  method . 

Among  the  molecular  orbitals,  the  la  and  2c  states  are 
practically  F Is  and  Li  Is  core  electron  states  respectively. 
Therefore,  their  wave  functions  are  calculated  only  in  their 
corresponding  atomic  regions  which  for  this  purpose  extend 
beyond  their  spherical  boundaries. 

The  total  energy  of  the  cluster  converges  more  slowly 
than  the  one-electron  energies,  as  the  self-consistent  iter- 
ations proceed.  The  total  energy  of  a LiF  molecule  is  about 
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1 

200  Ry,  and  it  has  to  be  conputsd  accurately  to  the  order  of 

0.001  Ry  in  order  to  obtain  an  accurate  picture  of  the  change 

of  the  molecular  total  energy  as  a function  of  internuclear 

distance  near  its  equilibrium  value.  Such  an  accuracy  could 

be  achieved  by  performing  enough  iterations  which  runs  about 

18  19 

15  with  a modified  Pratt  improvement  scheme  ' which  were 
incorporated  in  the  original  program  for  the  present  study. 

At  the  end  of  the  iteration  procedure  the  potential  conver- 
gence factor,  which  is  the  largest  ratio  of  the  difference 

between  the  starting  and  resultant  potentials  to  the  starting 

-3  -4 

potential,  becomes  of  the  order  of  10  "v  10  . The  normali- 

zation of  the  partial-wave  components  for  electron  states  was 
also  made. 

B.  Results 

Table  II  shows  the  results  of  the  non-SC?  (atomic)  and 

SCF  calculations  for  the  LiF  molecular  cluster  by  the  Xa 

multiple-scattaring  method.  The  non-SCF  calculation  was  done 

with  the  origi.nal  ionic  radii  for  the  atomic  spheres  and  the 

starting  potential  calculated  from  the  charge  densities  of 

neutral  atoms.  Also  shown  in  the  table  is  the  result  of 

20 

LCAO-MO  calculation  by  Ransil. 

As  self-consistent  iterations  proceed,  the  electronic 
charges  tend  to  be  more  uniformly  distributed  throughout  the 
) cluster  regions.  Consequently  the  magnitudes  of  all  the 

statistical  energies  decrease,  resulting  in  a lower  total  energy 
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(more  negative)  and  a larger  cohesive  energy.  At  the  same 
time  the  constant  potential  in  the  interatomic  region  becomes 
higher  and  so  do  all  the  one-electron  energies.  These  are 
the  general  tendencies  of  SC?  iterations  for  any  cluster  model 
calculations;  In  general,  the  Xa  one-electron  energies  are 
higher  than  the  corresponding  Hartree-Fock (KF)  one-electron 
energies  for  either  molecular  clusters  or  atoms.  There  is  a 

. i 

basic  difference  in  the  meaning  of  one-electron  energies  be-  ■ 

7 j 

tween  the  Xa  and  KF  methods. 

We  first  study  the  molecule  with  its  internuclear  distance  | 

at  experimental  equilibrium  value.  Except  the  Ic  and  2a  core  ^ 

states,  the  only  bound  electron  states -of  the  molecule  are  3a, 

4a,  and  Irr  states  which  are  mostly  made  of  F 2s/  F 2p  , and 

3 

F 2p^-type  atomic  wave  functions  respectively,  with  small  con- 
tributions from  s and  p functions  of  the  Li  atom.  Their  wave 
functions  are  shown  in  Figs.  2a,  2b,  and  2c;  and  their  normalized 
partial-wave  components  and  electron  charge  distributions  are 
listed  in  Table  III.  The  3a  and  4a  wave  functions  show  nodes 
in  the  Li  atomic  center  corresponding  to  the  node  of  the  Li  2s 
atomic  wave  function.  This  apparent  antibonding  character  of 
the  4a  state  causes  the  rising  of  the  4a  level  above  the  la 
level,  thus  reversi.ng  the  order  of  the  corresponding  LCAO-MO 
result,  as  shown  in  Table  II. 


Next  the  SCF  calculations  were  pursued  for  homogeneously 
compressed  and  expanded  clusters  to  obtain  the  theoretical 
equilibrium  properties  of  the  molecule.  The  results  of  the 
calculations  at  several  internuclear  distances  are  shown 
in  Table  IV.  As  the  cluster  expands,  the  constant  potential 
in  the  interatomic  region  and  the  magnitudes  of  all  the 
statistical  energies  obviously  decrease.  Since  the  atoms  are 
more  separated  at  the  saune  time,  each  electron  energy  level 
of  the  cluster  approaches  the  corresponding  electron  energy 
level  of  a free  ion,  and  the  ionicity  of  the  molecule  slightly 
increases  over  the  range  shown. 

The  ionicity  of  each  atomic  center  was  calculated  as 
follows.  At  the  end  of  an  iteration  we  get  electronic  charge 
distribution  for  each  valence  electron  state  over  the  various 
regions  of  the  cluster.  Since  the  electronic  charges  in  the 
interatomic  and  extramolecular  regions  do  not  belong  to  any 
specific  atom,  we  reassign  them  to  the  atomic  centers  for 
each  valence  electron  state  separately,  in  proportion  to  the 
valence  electronic  charges  which  already  exists  in  the  atomic 
ceneters  as  a result  of  the  calculation.  Then,  the  sum  of 
all  the  charges  in  an  atomic  center  gives  the  ionicity  of  the 
atom  in  the  cluster.  This  redistributions  of  the  electronic 
charges  is  done  only  to  calculate  the  ionicities. 

The  virial  coefficient  is  defined  here  as  the  ratio  of 
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the  total  idnetic  energy  to  the  magnitude  of  the  total  energy. 

For  all  the  cluster  configurations  treated,  the  virial  co- 
efficients are  very  close  to  unity  because  the  overlapping- 
sphere  scheme  was  used  such  that  the  virial  theorem  can  be 
closely  satisfied.  The  coefficient  becomes  smaller  as  the 
atomic  centers  are  overlapped  more  for  a given  internuclear 
distance,  however,  more  overlapping  than  used  here  is  not 
pursued  for  the  reason  stated  earlier.  The  virial  coefficient 
becomes  closer  to  unity  for  more  expanded  cluster,  since 
the  system  then  is  more  like  a collection  of  free  ions  for  which 
the  virial  theorem  is  better  satisfied. 

The  total  energy  of  the  LiF  molecule  is  plotted  against 
the  internuclear  distance  in  Fig.  3.  The  figure  shows  that 
there  is  a stable  equilibrium  state  over  the  range.  We  now 
use  a simple  parabolic  formula  for  the  total  energy  as  a 
function  of  internuclear  distance, 

E(R)  = J / (1) 

and  fit  the  total  energies  at  the  latter  three  internuclear 

distance  of  R , 1,02  R , and  10.5  R , where  R^  is  the 
e e e e 

experimental  equilibrium  internuclear  distance  of  the  molecule. 

The  E is  the  total  energy  of  the  molecular  cluster  at  the 
o 

theoretical  equilibrium  internuclear  distance  R^,  and  the 
constant  k then  can  be  interpreted  as  a spring  constant  be- 
tween two  atomic  centers  of  the  molecule. 
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Assuming  the  harmonic  motion  for  the  vibration  of  the 
molecule  at  the  equilibrium  configuration,  the  angular  frequency 
0)^  of  the  vibration  is  given  by 

0)^  * , (2) 

where  u is  the  reduced  mass  of  the  molecule. 

The  properties  of  the  LiF  molecule  thus  calculated  from 
the  model  are  listed  in  Table  V along  with  corresponding  ex- 

t 

perimental  data.  The  theoretical  equilibrium  internuclear  dis- 
tance is  only  about  1.3%  larger  than  the  experimental  value, 
and  the  calculated  vibrational  quantum  energy  is  also  in  very 
good  agreement  with  the  experimental  value. 

The  cohesive  energy  was  obtained  as  the  difference  be- 
tween the  total  energy  of  the  molecular  cluster  and  the  sum 
of  the  Xci  total  energies  of  a Li  and  an  F atoms.  The  exchange 
a parameters  for  atom.s  were  computer  by  Schwarz^"^  such  that 
the  Xa  total  energy  of  an  atom  is  equal  to  the  restricted 

Hartree-Fock  total  energy  of  the  atom  which  was  listed  i.n  a 
21 

report  by  .Mann  for  all  atoms.  In  Table  V the  calculated  total 
energy  is  smaller  in  magnitude  than  the  experimentally  observed 
total  energy  while  the  calculated  cohesive  energy  is  larger 
than  the  experimental  cohesive  energy.  This  is  due  to  the 
difference  in  the  atomic  total  energies  between  the  theoretical 
and  exoerimental  situations.  The  cohesive  energies  of  the 
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of  the  molecule  in  both  cases  are  reasonably  close,  however. 

The  theoretical  electric  dipole  moment  was  obtained  as 
the  product  of  the  calculated  values  of  the  ionicity  and  the 
equilibrium  internuclear  distance.  It  is  in  very  good  agree- 
ment with  the  experiment. 

There  is  one  renark  to  be  made  concerning  the  accuracy 
of  the  calculations.  As  we  see  from  Fig.  3,  there  is  a rather 
large  margin  of  uncertainty  for  the  total  energy  of  the  mole- 
cular cluster.  (Also  see  previous  subsection  A.)  More 
accuracy  can  not  be  attained  since  the  relative  error  of  the 
total  energy  at  the  last  stages  cf  iterations  is  of  the 
order  of  10  ^ 10  ^ which  is  close  to  the  minimum  error 

possible.  The  uncertainty  affects  the  vibrational  quantum 
energy  most  seriously,  which  has  to  be  interpreted  with  the 
margin  of  error  of  about  +20%. 
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III.  LiF  CRYSTAL 

The  LiF  in  a crystalline  form  is  an  insulator  with  a 

22 

direct  band  gap  of  1.000  + 0.004  Ry  . It  is  of  a NaCl 

structure  and  its  constituent  atoms  are  the  lightest  among 

alkali  halides.  Consequently  it  has  the  smallest  equilibrium 

nearest-neighbor  distance  among  the  family  of  crystals,  which 

* 23 

is  3.998  a.u.  (2.01  A)  being  exactly  the  sum  of  the  Zachariasen ' s 
ionic  radii  for  Li  and  F as  given  in  Table  I.  The  nearest- 
neighbor  distance  for  LiF  has  been  increased  by  about  25%  as 
the  molecules  form  a crystal  primarily  due  to  the  isotropic 
electrostatic  interactions  of  ions  in  the  crystalline  environ- 
ment. In  the  crystalline  form,  the  LiF  will  still  maintain  a 
high  ionicity. 

A.  Cluster  Model 

If  we  choose  the  smallest  possible  cluster  model  for  the 

LiF  crystal  while  retaining  the  full  point  group  symmetry  of 

the  crystal,  it  is  (LiFg)”^  or  (FLig)^^  necessarily  including 
24 

a Watson  sphere  in  order  to  make  the  whole  system  neutral. 
However,  they  are  too  much  unbalanced  in  composition  of  atoms 
to  simulate  the  LiF  crystal  properly.  These  cluster  models 
are  also  unsuitable  for  the  F-center  problem  of  the  crystal 
which  will  be  considered  in  Section  IV. 

If  we  add  two  more  atomic  layers  of  atoms  to  (FLig)"*"^, 
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for  example,  the  cluster  structure 


becomes  (FLig 


which  is  almost  neutral  but  an  extra  electron  still  requires 
the  use  of  a Watson  sphere.  This  could  be  a good  cluster 
model  for  the  LiF  crystal  and  also  for  the  F-center  problem. 
However,  it  has  two  inequivalent  atomic  sites  for  each  of  Li 
and  F atoms  unlike  the  crystal.  It  is  also  rather  large  for 
the  present  purpose  of  simulating  the  crystal  by  a conveniently 
small  cluster  model.  . . 

The  smallest  neutral  cluster  model  for  the  LiF  crystal 
larger  than  a mere  molecule  is  a cluster  iri  which  the  atomic 

centers  are  placed  tetrahedrally  at  the  corners  of  a cube  for 
each  atomic  species.  The  next  layer  of  atoms  contains  24  atoms 
altogther  which,  if  included,  would  make  the  cluster  quite 
large.  In  the  Li^F^  cluster  model,  all  the  atoms  of  each  species 
are  equivalent  as  they  are  in  LiF  crystal.  The  compositions 
are  also  the  same  in  both  the  cluster  and  crystal. 

In  general  a balanced  cluster  for  a compound  is  less 
symmetric  than  the  crystal;  the  Li^F^  cluster  has  the  T^^  point 
group  symmetry  instead  of  the  point  group  symmetry  of  the 
LiF  crystal.  However,  the  reduction  of  symmetry  will  not  by 
itself  present  a serious  problem  because  we  can  always  use 
the  compatibility  relations  for  representations  between  a 
symmetry  group  and  its  subgroup  to  identify  electron  states. 
Furthermore ,. in  any  cluster  model  regardless  whether  it  has 
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a full  point  group  symnetry  of  the  crystal  or  not,  only  the 
states  of  the  highest  syrmetry  representations  can  be  related 
to  the  states  of  the  crystal  if  it  is  ever  possible. 

Since  the  cluster  is  neutral,  the  electrostatic  inter- 
action between  the  cluster  and  the  remainder  of  the  crystal 
is  kept  at  a minimum  for  the  given  size  of  the  cluster. 

The  results  of  the  SCF  Xa  claculations  for  the  Li.F. 

4 4 


cluster  with  various  nearest-neighbor  distance  within  ^20% 
of  the  experimental  equilibrium  value  revealed  no  equilibrium 
state  for  the  cluster  model.  Only  when  an  empty  spherical 
center  (call  0 center)  was  added  at  the  geometric  center  of 
the  cluster  could  an  equilibr i'u.m  state  be  reached  at  a 
nearest-neighbor  distance  not  far  from  the  experimental  equili- 
brium value.  The  empty  0 center  is  treated  just  like  an  atomic 
center  in  the  multiple-scattering  calculations,  and  its  in- 
clusion in  the  cluster  reduces  the  interatomic  region  where 
wave  functions  are  treated  less  accurately  than  in  the  atomic 
regions.  Also  the  bonding  behaviors  between  the  atomic  centers 
will  be  described  better  by  linking  the  scattered  waves  from 
the  atomic  centers  more  closely,  especially  because  the  0 center 
takes  up  the  central  position  among  the  atomic  centers.  The 
model  will  be  called  OLi^F^  cluster  model,  and  we  use  it  to 
investigate  the  properties  of  the  LiF  crystal  and  later  study 
the  F-center  problem.  The  cluster  model  is  shown  in  Fig.  4. 
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B. 


Calculational  Considerations 


In  the  SCF  Xa  multiple-scattering  calculations  for  the 
OLi^F^  cluster,  the  same  exchange  parameter  a • s of  the  LiF 
molecular  cluster  were  used  in  the  corresponding  regions  of 
*the  cluster.  The  a in  the  0 center  is  the  same  a used  in  the 
interatomic  and  extramolecular  regions. 

The  symmetrized  basis  functions  for  the  valence  electron 
states  are  made  of  the  s and  p functions  of  Li  and  F together 
with  the  s,  p,  d,  and  f functions  for  the  0 center  and  the 
extraumolecular  region. 

Since  the  size  of  the  atomic  centers  are  more  restricted 
by  the  presence  of  larger  number  of  neighbors  in  the  OLi^^F^ 
cluster  than  in  the  molecular  cluster,  the  radii  of  the  atomic 
centers  were  calculated  more  rigorously  as  follows  for  the 
cluster  model  with  experimental  equilibrium  nearest-neighbor 
distance.  First,  the  touching  radii  of  the  Li  and  F atomic 
centers  were  determined  such  that  the  spherically  averaged 
potentials  in  any  neighboring  atomic  centers  are  equal  at 
the  contact  point  of  the  atomic  spheres.  These  touching  radii 
remained  essentially  constant  after  only  a few  iterations;  they 
are  1.675  a.u.  (0.866  A)  for  Li  and  2.123  a.u.  (1.123  A)  for 
F atomic  center,  being  intermediate  between  the  atomic  and 
ionic  radii.  The  touching  radius  of  the  0 center  is  1.166  a.u. 
(0.619  A)  which  is  about  50%  larger  than  the  corresponding 
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radius  when  the  ionic  radii  are  used  for  the  Li  and  F atomic 
centers.  In  the  subsequent  interations , then,  the  touching 
radii  are  linearly  expanded  so  that  the  virial  theorem  is 
satisfied  for  the  cluster.  In  this  process  the  potentials  of 
the  neighboring  atomic  centers  at  the  touching  radii  still 
remain  practically  equal.  This  procedure  results  in  an  over- 
lapping scale  factor,  1.14,  and  it  is  used  for  the  clusters 
with  different  nearest-neighbor  distances  as  well. 

The  Madelung  corrections  to  the  atomic  potentials  are 
included  to  account  for  the  neglected  electrostatic  inter- 
action between  the  ions  in  the  cluster  and  the  ions  in  the  re- 
mainder of  the  crystal  which  were  excluded  from  the  cluster 
model.  The  Madelung  correction  term  to  the  potential  in 
the  th  atomic  center  is  calculated  by 

1 q- 

. ^rt^l  ion 


where  R is  the  nearest-neighbor  distance,  and  is 

ionicity  of  the  £ th  atomic  center  which  is  calculated  i.n 

the  same  way  as  described  in  the  Subsection  3 of  the  previous 

Section  II.  That  is,  we  assign  the  ionicity  of  the  atomic 

centers  of  the  cluster  to  the  ions  in  the  remainder  of  the 

c L 

crystal.  The  .Madelung  constant  of  the  cluster  a „ was  obtained 
by  subtracting  the  contributions  of  the  ions  within  the  cluster 
from  the  .Madelung  constant  of  the  crystal;  this  turns  out  to 


be  a 0.2915  which  is  significantly  reduced  from  the  Madelung 
constant  of  the  crystal,  1.7476.  In  the  present  treatment 
the  Madelung  correction  terms  to  the  potentials  in  the  inter- 
atomic and  extramolecular  regions  as  well  as  the  central  0 
region  vanish. 

C . Results  and  Discussions 

Table  VI  lists  the  results  of  the  SCF  Xx  multiple-scattering 

calculations  for  the  OLi.F,  cluster  with  various  nearest- 

4 4 

neighbor  distances.  They  show  the  same  tendencies  as  the 
cluster  exoands  as  we  discussed  earlier  with  the  LiF  molecular 


cluster . 

The  total  energy  of  the  cluster  is  plotted  against  lattice 
constant  (which  is  twice  the  nearest-neighbor  distance  in  LiF 
crystal)  in  Fig.  5.  The  equilibrium  state  of  the  OLi^F^  cluster 
is  attained  whe.n  the  lattice  constant  is  7.4%  larger  than  the 
experimental  equilibrium  value.  This  result  is  satisfactory 
considering  the  size  of  the  cluster  used.  At  equilibrium 
the  total  energy  of  the  cluster  is  -214.322  Ry  which,  compared 
with  the  sum  of  the  total  energies  of  free  atoms,  gives 
cohesive  energy  per  molecule,  0.638  Ry.  This  is  in  an  excellent 
agreement  with  an  experimental  value,  0.644  Ry.  The  experi- 
mental cohesive  energy  with  respect  to  free  atoms  was  deduced 

25 

from  the  cohesive  energy  with  respect  to  free  ions,  -0.737  Ry, 
with  corrections  made  for  the  ionization  energy  of  a Li  atom. 
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0.396  Ry/  and  the  electron  affinity  of  an  F atom,  -0.253  Ry. 

Out  of  the  calculated  cohesive  energy,  about  0.1  Ry  was  con- 
tributed from  the  Madelung  correction. 

Now  in  order  to  calculate  the  bulk  modulus  we  again  assume 
a parabolic  formula  of  Eq.  (1)  for  the  total  energy  of  the 
cluster  near  equilibrium,  and  obtain  a spring  constant  k by 
fitting  the  total  energies  at  the  three  nearest-neighbor 
distances  R = 1.05  Re,  1.074  Re,  and  1.10  Re,  where  Re  is  the 
experiemental  equilibrium  nearest-neighbor  distance.  Then 
the  equilibrium  bulk  modulus  3^  is  given  by 
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where  V is  the  volume  for  four  molecules  in  the  crystal,  i.e. 

V = 8R^,  and  subscript  o indicates  that  the  parameters  are 

taken  from  the  cluster  at  equilibrium.  Thus  calculated,  the 

11  , 

equilibrium  bulk  modulus  6.66  x 10  dyne/cm“  is  also  in  very 

11  2 

good  agreement  with  an  experimental  value  6.93  x 10  dyne/cm 

which  was  obtained  from  the  elastic  stiffness  constants  of 

2 8 

Brisco  and  Squire. 

The  ionicity  of  the  LiF  crystal  through  the  OLi^F^ 
cluster  model  is  0.399  at  equilibrium  in  units  of  electronic 
charge.  The  reduction  of  the  ionicity  from  0.92  for  the  LiF 
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molecule  is  seen  primarily  due  to  the  increased  nearest- 

neighbor  distance.  The  quilibrium  properties  of  the  LiF 

crystal  obtained  by  the  SCF  Xa  multiple-scatteri.ng  calculations 

for  the  OLi.F,  cluster  model  are  summarized  in  Table  VII. 

4 4 

The  one-electron  energy  levels  of  the  equilibrium  OLi^F^ 
cluster  are  shown  in  Fig.  6(a).  The  levels  for  the  core 


electrons  and  mostly 


F 2s  atomic  electrons  are  not  shown. 


In  the  figure  the  set  of  valence  electron  levels  corresponds 
to  the  top  valence  band  of  the  LiF  crystal  while  the  excited 
levels  correspo.nd  to  the  conduction  band  levels  of  the  crystal. 

In  fact,  the  symmetries  of  the  3t2  and  3a^  states  conform  to 
the  symmetries  of  the  top  electron  state  of  the  valence 

band  and  the  bottom  electron  state  of  the  conduction  band 
of  the  crystal  respectively  through  the  com.patibility  relations 
between  the  T^  and  0^  point  groups. 

The  partial-wave  coefficients  and  the  charge  distributions 
of  the  3t2,  3aj^,  and  4t2  states  are  described  in  Table  VIII, 
and  the  contour  maps  of  their  wave  functions  are  shown  in  Figs. 

7a,  7b,  and  7c.  The  3t2  state  is  mostly  composed  of  an  F 2p 
atomic  basis  state,  but  with  a small  contribution  from  a Li  2s 
atomic  basis  state  it  gives  a slight  bonding  between  F and  Li 
atoms.  The  Sa^^  state  is  in  a large  part  an  extramolecular  s state 
having  a little  more  than  half  of  its  electronic  charge  pro- 
bability in  the  extramolecular  region.  The  higher  4t2  state 


3 


21 


I 


) 

i 

i 


has  even  more  of  its  electronic  charge  probability  in  the 
extramolecular  region,  primarily  as  an  extraumolecular  p state. 

The  excited  states,  la^^  and  ^t^,  have  more  or  less  uniform 
distribution  of  the  probability  of  finding  electrons  through- 
out the  regions  of  the  cluster.  This  is  characteristic  of 
the  conduction  band  electron  states  of  a crystal. 

The  bandgap  derived  from  the  cluster  model  is  0.54  Ry.  It 
is  quite  small  com.pared  to  the  bandgap  of  the  LiF  crystal  (about  1 Ry)  . 
This  discrepancy  can  be  explained  as  follows.  In  the  crystal 
the  conduction  band  levels  are  pushed  up  by  the  repulsive 
interactions  from  the  electrons  which  has  already  filled  the 
valence  band  levels  com.pletely  and  are  located  mostly  in  the 
atomic  regions.  I.n  the  cluster  model,  however,  the  electrons 
of  the  excited  states  largely  escape  from  the  atomic  regions 
to  the  extramolecular  region  where  little  valence  electronic 
charge  resides.  Consequently  the  excited  state  electrons  can 
be  more  comfortably  accommodated  in  the  cluster  than  in  the 
crystal,  thus  lowering  their  energy  levels  below  where  they 
would  be  in  the  crystal.  The  width  of  the  valence  band  of  the 
cluster  is  also  significantly  more  narrow  than  the  width  of 
the  corresponding  valence  band  of  the  crystal  because  the 
interactions  a.mong  the  valence  electrons  in  the  cluster  are 
limited  by  the  presence  of  the  empty  extramolecular  region, 
that  is,  by  omitting  some  nearest-neighbor  atoms. 


) 


Now  suppose 
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that  all  the  energy  levels  of  the  valence  electrons,  including 

the  excited  levels,  were  contracted  by  the  sane  degree  in 

the  cluster  model  by  the  presence  of  the  extramolecular  region. 

Then,  in  order  to  compare  our  results  with  energy  band  of  the 

LiF  crystal,  we  linearly  expand  the  energy  levels  of  the 

cluster  in  the  ground  state  configuration  such  that  its  band- 

gap  is  equal  to  1 Ry,  the  bandgap  of  the  crystal,  while 

aligning  the  tops  of  the  valence  bands.  The  linear  scale 

factor  is  1.35.  The  width  of  the  resultant  valence  band  of 

the  cluster  is  0.137  Ry,  now  quite  close  to  the  valence  band 

width,  about  0.125  Ry,  of  Brener's  energy  bands  for  the  LiF 
29 

crystal.  The  linearly  expanded  energy  levels  of  the  cluster 
are  shown  in  Fig.  8(a),  compared  with  the  Brener's  energy 
bands  in  Fig.  8 (b) . 

In  order  to  find  an  exciton  level  in  the  present  model 
we  consider  a transition-state  cluster,  using  Slater's 
transition-state  concept,"^  which  has  half  an  electron  trans- 
ferred from  the  top  valence  state, 3 , to  the  lowest  excited 
state,  3a^.  The  resultant  energy  levels  of  this  cluster  is 
shown  in  Fig.  6(b).  Here  the  energy  level  of  the  3a state 
is  lowered  by  0.044  Ry  from  the  corresponding  level  in  the 
ground-state  configuration.  It  seems  that  the  la^^  stats  in 
the  transistion-state  configuration  represents  the  first  exciton 
state,  while  the  la^^  state  in  the  ground-state  configuration 
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represents  the  bottom  of  the  conduction  band  of  the  crystal. 

* 22 
Experimentally  it  was  shown  by  Roessler  and  Walker  that  the 

first  exciton  level  in  the  fundamental  absorption  of  the  LiF 

crystal  is  at  0.927  Ry  (with  the  half-width  of  0.024  Ry) , 

i.e.  0.073  Ry  below  the  band-to-band  transition.  When  the 

energy  levels  of  the  transition-state  cluster  are  also  expanded 

by  the  sa.me  linear  scale  factor,  the  exciton  level  is  placed 

0.081  Ry  below  the  bottom  of  the  conduction  band.  This  is 

now  in  quite  good  agreement  with  the  experiment. 
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IV.  F-CENTER  IN  THE  LiF  CRYSTAL 

The  F-centers  in  alkali  halides  have  been  the  subjects  of 
many  extensive  investigations  largely  due  to  the  significance 
in  optical  properties,  and  they  are  reviewed  by  Markham^® 
and  Fowler. 

On  the  theoretical  side  of  the  approach  to  the  F-center 
problem,  there  have  been  calculations  using  a modified 
Hartree-Fock  equation.  One  is  t;he  point-ion  approximation 
based  on  neglecting  the  orthogonalization  of  the  defect 
electron  states  to  other  electron  states.  The  results  ob- 
tained by  a variational  technique  give  absorption  energies 
of  the  F-centers  in  alkal  halides  to  within  10  'v  15%  of 

the  experimental  values.  Some  improvements  have  been  made 

32 

by  the  exte.nded-ion  approximation  and  the  pseudo-potential 
approximation.^^  More  crude  semi-continuum  model  is  also 
satisfactory  for  the  F-center  absorption  problem.  In  fact, 
just  about  any  model  involving  a potential  well  gives  rea- 
sonably good  results  for  the  absorption  energies. 

There  is  a recent  first-principle  LCAO  calculation  by 
Chaney  and  Lin^^  for  the  electronic  structure  of  the  F-center 
in  the  LiF  crystal.  In  this  calculation  atomic  orbitals  are 
approximate  as  being  Gaussian-type  to  facilitate  the  multi- 
center integrals  associated  with  the  Hamiltonian  matrix  elements. 
The  multiple-scattering  method  provides  another  first- 
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principle  approach  to  the  problem,  and  it  is  simpler  and  more 
straight-forward  while  keeping  the  orthogonality  of  the 
electron  states  in  the  system.  There  is  no  need  to  calculate 
difficult  multicenter  integrals,  therefore,  the  SCF  cal- 
culations can  be  performed  rapidly.  In  this  case  approximation 
is  made  in  modeling  a cluster  suitable  for  the  crystal  con- 
tinaing  the  vacancy. 

The  SCF  Xa  multiple-scattering  method  can  be  applied  to 
a large  cluster  for  complete  analysis  of  the  F-center  problem. 
However,  we  do  not  intend  to  be  complete.  Rather  we  will 
start  from  the  small  cluster  model  for  the  perfect  LiF  crystal, 
and  show  how  much  information  we  can  obtain  with  such  a min- 
imal size  of  the  cluster  containing  the  F-vacancy. 

A.  F-Vacancy  Model  and  Calculaticnal  Considerations 

In  formulating  a cluster  model  for  the  F-vacancy  in  the 
LiF  crystal,  we  will  be  consistent  with  the  OLi^F^  cluster 
modeled  for  the  perfect  crystal  so  that  we  can  also  compare 
the  results  of  both  cases  directly.  Therefore,  we  remove 
the  F atom  from  an  F atomic  center  in  the  OLi^F^  cluster,  and 
the  resultant  empty  F atomic  center  (now  called  VacF  center) 
is  just  like  the  0 center.  This  F-vacancy  cluster  model 
becomes  a cluster  of  O-VacF-ALi-Li^-F^  (called  VacF  cluster 
for  convenience) . It  has  a unique  symmetry  axis  which  connects 
•the  geometric  centers  of  the  VacF  center  and  the  0 center,  and 
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the  ALi  atomic  center  is  a Li  atomic  enter  on  the  axis, 
thereby  being  inequivalent  to  the  other  three  Li  atomic  centers, 
which  are  equivalent  among  themselves.  The  remaining  three 
F atomic  centers  are  also  equivalent.  The  symmetry  of  the 
VacF  cluster  belongs  to  the  point  group,  a subgroup  of 

the  T^  point  group  for  the  OLi^F^  cluster. 

In  the  SCF  Xa  multiple-scattering  calculations  for  the 
VacF  cluster,  the  same  .Madelung  correction  as  was  used  for 
the  OLi^F^j  cluster  is  kept  with  the  same  ionicity.  The  ex- 
change a parameter  for  the  VacF  center  is  the  same  a as  we 
adopted  for  the  0 center  as  well  as  for  the  interatomic  and 
extramolecular  regions  of  the  OLi^F^  cluster.  Otherwise, 
the  exchange  a parameters  and  geometrical  factors  of  the  OLi^F^ 
cluster  in  its  equilibrium  state  are  retained. 

The  defect  electron  states  due  to  the  F-vacancy  will  be 
described  in  the  extramolecular  regions  where  neighboring 
atoms  are  missing  as  well  as  in  the  atomic  centers  and  the 
interatomic  region.  In  that  sense  the  VacF  cluster  model  is 
not  quite  ideal,  however,  it  will  be  the  smallest  possible 
neutral  cluster  model  for  the  F-vacancy,  and  in  the  present 
framework  the  defect  electron  states  are  equally  well  treated 
as  the  valence  electron  states  of  the  VacF  cluster  or  the 

OLi.F,  cluster. 

4 4 

The  electron  states  of  the  VacF  cluster  and  the  OLi^F^ 
cluster  can  be  correlated  through  tbe  use  of  the  compatibility 
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relations  between  the  representations  of  the  and 
point  groups  as  shown  in  Table  IX. 

B.  Results  and  Discussions 

The  VacF  cluster  is  an  open-shell  system  obviously  be- 
cause the  total  number  of  electrons  are  odd,  being  39  in  all. 
Therefore,  a spin-polarized  scheme  is  in  order  to  calculate 
electron  states  for  the  cluster.  Calculations  with  spin- 
restricted  scheme  for  the  VacF  cluster  and  (VacF)"*”^  ionic 
cluster  (the  one  missing  electron  from  the  F-center  is  now 
frozen  on  a Watson  sphere,  which  is  geometrically  the  same  as 
the  outersphere,  to  make  the  total  system  neutral)  are  also 
done  for  a comparative  study.  The  results  for  the  various 
clusters  are  described  in  Table  X,  and  their  electron  energy 
level  structures  are  shown  in  Fig.  9 compared  with  that  of 

the  OLi.F.  cluster. 

4 4 

As  an  F atom  is  removed  from  the  OLi^F^  cluster,  the  top 
valence  band  loses  six  electron  states  counting  degeneracy 
which  accounts  for  a missing  six-fold  degenerate  F 2p  atomic 
level.  At  the  same  time  only  five  valence  electrons  are 
lost,  therefore,  one  electron  is  left  over  after  the  valence 
band  is  completely  filled.  This  extra  electron,  which  is 
associated  with  the  F-center,  is  consequently  accomodated  by 
the  next-higher  4aj^  state  in  the  VacF  cluster.  It  is  largely 
contained  in  the  VacF  center  in  the  VacF  cluster  as  one  can 


see  by  comparing  the  ionicities  of  the  VacF  center  in  the 
VacF  cluster  and  the  (VacF)'*’^  cluster. 


In  the  transition  of  the  OLi.F.  cluster  to  the  (VacF)"*"^ 

4 4 

cluster  by  losing  an  f”  ion,  the  width  of  the  valence  band 
decreases  from  0.074  Ry  to  0.065  Ry  since  six  valence  electron 
states  are  lost.  At  the  same  time  the  lost  six  electrons  from 
cluster  account  for  the  lowering  of  the  valence  and  conduction 
band  levels  by  decreasing  electron-electron  interactions  for 
the  remaini.ng  electrons.  Particularly,  the  4a^  level  is 
lowered  deep  into  the  bandgap  from  the  conduction  band  by  con- 
centrating its  wave  function  in  the  VacF  center  which  would 
acconmodate  excited  electrons  .more  comforatably  than  a.ny  other 
region  of  the  cluster. 

The  introduction  of  an  electron  into  the  (VacF)'*’^  cluster, 
thus  freeing  the  electron  on  the  Watson  sphere  to  make  the 
system  the  VacF  cluster,  raises  the  valence  band  levels  by 
increasing  electron-electron  interaction.  However,  the  con- 
duction band  states  are  lowered  a little  because  the  new  electron 
moves  towards  the  inside  of  the  cluster  away  from  the  Watson 
sphere  around  which  the  wave  functions  of  the  conduction  band 
states  are  large. 

Next,  spin-polarization  of  the  VacF  cluster  lowers  all 
the  spin-up  electron  levels  (the  4aj^  defect  electron's  spin 
is  arbitrarily  assu-med  to  be  up)  and  raises  all  the  spin- 


down  electron  levels,  though  changes  in  the  valence  band. are 
negligible.  The  energy  level  structure  of  the  spin-down 
electrons  of  the  VacF  cluster  is  quite  similar  to  that  of 
the  OLi^F^  cluster.  For  all  the  VacF  cluster  configurations 
the  width  of  the  balance  band  remains  fairly  constant. 

The  wave  functions  of  the  top  valence  band  state,  4e, 
and  the  excited  states,  4aj^,  5a^,  and  5e,  for  both  spins  are 
shown  in  Figs.  lOa^  through  10d2  by  their  contour  map;  and 
their  norm.alized  partial-wave  components  and  electronic  charge 
distributions  are  listed  in  Tables  XIa  and  Xlb. 

The  4e  and  4e  states  are  almost  identical  to  each 
other  and  similar  to  the  top  valence  state,  3t_,  of  the  OLi.F. 
cluster,  except  in  the  centers  along  the  symmetry  axis  including 
the  VacF  center.  In  terms  of  a spin-restricted  perturbation, 
the  3tT  state  of  the  OLi .F.  cluster  is  SDlit  into  an  a,  and  an 
e states  as  its  symmetry  is  reduced  from  the  T^  to  point 

group  with  the  introduction  of  the  VacF  center.  The  e state 
of  the  F-vacancy  retains  most  of  the  characters  of  the  3t2 
state  of  the  OLi^F^  cluster  except  along  the  symmetry  axis. 
However,  the  a^  state  which  had  its  wave  function  concentrated 
along  the  axis  previously  in  the  OLi^F^  cluster  is  now  so 
much  perturbed  by  the  F-vacancy  that  its  energy  level  is  driven 
up  enormously  (not  shown  in  Fig.  9c)  . This  a^^  state  together 
with  another  e'  state  from  the  valence  band  accounts  for  the 
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six  electron  states  lost  from  the  band  as  we  discussed 
earlier. 

The  wave  function  of  the  4aj^^  state  is  largely  centered 
around  the  F-vacancy  region  as  we  see  from  Fig.  lOb^,  thereby 
simulating  the  F-center  defect  electron  state  of  the  LiF 
crystal.  Its  electronic  charge  is  fairly  uniformly  distri- 
buted in  the  VacF  center  with  the  total  of  about  half  an 
electron.  The  wave  function  is  mostly  s-type  with  most  of 
the  contributiions  coming  from  the  VacF  s and  extramolecular 
s basis  functions.  The  4a^  states  of  both  spins  are  the  per- 
turbed states  of  the  Sa.  state  of  the  OLi.F,  cluster  due  to  the 

1 4 4 

F-vacancy . 

+ 4 

The  diff  erence  - between  the  4aj^  and  4aj^  states  is  seen 
to  result  mainly  from  the  difference  in  potentials  in  the  F- 
vacancy  region  between  spin-up  and  spin-down  electrons,  as 
shown  in  Fig.  11.  The  potential  for  a spin-up  electron  is 
rather  low  and  approximately  constant  at  about  -0.65  Ry.  On 
the  other  hand  the  potential  for  a spin-down  electron  becomes 


The  4t2  state  of  the  OLi^F^  cluster  is  now  split  to  the 
5e  and  5a^  states  when  the  F-vacancy  is  introduced.  Let  us 
first  look  at  the  5aj_^  and  5a states.  They  are  similar  to 
each  other  in  many  aspects.  But,  while  the  primary  extramole- 
cular  p character  of  the  original  4t2  state  of  the  OLi^F^ 

L 

cluster  is  well  preserved  by  the  Sa^^  state  of  the  VacF  cluster, 
a significa.nt  amount  of  s function  sets  i.n  to  distort  such 
p character  in  the  extramolecular  region  for  the  5aj^  state. 

The  extramolecular  s component  is  an  extension  of  the  s-type 
fucntion  of  the  5a^^  state  in  the  VacF  center  which  also  ex- 
tends to  the  nearest-neighbor  Li  atomic  centers  of  the  cluster. 
Such  character  of  the  Sa^^  wave  function  proves  that  the  extra- 
molecular  region  can  approximately  simulate  missing  neighboring 
atomic  centers  of  the  F-vacancy  in  the  present  cluster  model. 

The  s-type  wave  functions  centered  around  the  F-vacancy  for  the 
Sa^^  states  are  less  localized  than  those  for  the  4aj^  states. 

The  5e  and  5e  states  are  quire  identical  to  each  other, 
and  mostly  made  of  extramolecular  p basis  functions  as  in  the 
case  for  the  original  4t2  state  of  the  OLi^F^  cluster.  The 
contour  map  of  the  4t2  wave  function  of  the  OLi^F^  cluster, 
which  is  shown  in  Fig.  7c,  has  p_-type  basis  fu.nction  in 
the  extramolecular  region.  We  can  always  combine  this 
wave  function  with  the  other  two  equivalent  P^-  and  p^- 
types  to  obtain  new  wave  functions  whose  extramolecular  p 
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functions  are  directed  differently,  for  example,  along  the  new 
symmetry  axis  of  the  VacF  cluster  as  its  5e  wave  functions. 
Similar  agruments  can  be  made  also  for  the  Sa^^  wave  functions. 
For  the  Sa^^  and  5e  states  also,  the  energy  differences  between 
the  spin-up  and  spin-down  electrons  are  attributed  primarily 
to  the  difference  in  potentials  in  the  F-vacancy  region. 

As  the  potential  for  a s?in-up  electron  suggests,  the 
F-center  proble.m  for  the  defect  electron  corresponds  quali- 
tatively to  a three-dimensional  potential-well  problem.  In 

34 

view  of  a simple  infinite  square  well  model  or  a hydrogenic 
model  for  the  F-centers  in  alkali  halides,  the  4a^  state  of 
■ the  VacF  cluster  corresponds  to  the  Is  state,  and  the  Sa^^^  and 
5e^  states  correspond  to  the  2s  and  2p  states  respectively.  In 
these  models  the  optical  absorption  of  the  F-center  is  the 
result  of  an  electron  excitation  from  the  Is  to  the  2p  state. 

A X 

The  excitation  energy  from  the  4aj^  to  the  5e  state  in 
the  VacF  cluster  model  is  0.178  Ry  as  the  energy  difference 
between  the  two  energy  levels.  The  transition  energy  calculated 
by  Slater's  transition-state  scheme  is  0.175  Ry,  not  much 
different  from  the  above  value.  The  difference  in  total  energies 
between  the  clusters  in  ground  and  excited  state  configurations 
gives  also  si.milar  figure  for  the  excitation  energy.  However, 
if  we  expand  the  energy  levels  of  ground-state  VacF  cluster 
linearly  by  the  same  scale  factor  of  1.85  as  we  did  for  the 
energy  levels  of  the  OLi^F^  cluster,  the  excitation  energy  be- 


comes  0.33  Ry,  close  to  an  experimental  result  which  indicates 
that  the  F-center  absorption  in  a LiF  crystal  is  peaked  at 
0.382  Ry  with  the  half-width  of  0.044  Ry.^^ 

In  the  transition  from  the  (VacF)'*’^  cluster  to  the  spin- 
restricted  VacF  cluster,  total  energy  increases  (in  magnitude) 
by  0.181  Ry.  Spin-polarization  in  the  VacF  cluster  makes  the 
total  energy  still  larger  (also  in  magnitude)  by  0.033  Ry  more. 

The  formation  energy  of  the  F-center  in  LiF  is  0.674  Ry 
in  the  cluster  model  as  the  difference  in  total  energies 
between  the  OLi_^F^  cluster  and  a system  of  the  VacF  cluster 
plus  a free  F atom.  However,  the  experimental  formation 
energy  is  given  with  respect  to  an  F atom  adsorbed  on  the 
surface  of  the  crystal;  it  is  0.172  ^^0.197  Ry.'^^  Because  such 
adsorption  energy  is  not  available,  a comparison  between  the 
cluster  model  calculation  and  experiment  for  the  formation 
energy  of  an  F-ce.nter  can  not  be  made  at 


present . 


V.  SUMMARY 


A LiF  moleculavas  first  examined  by  direct  application  of 
the  SCF  Xa  multiple-scattering  method,  and  the  results  for 
molecular  properties  were  shown  in  excellent  agreement  with 
experiments.  Then  we  applied  the  method  to  the  smallest- 
possible  neutral  cluster  model  of  OLi^F^ , consisting  of- only 
eight  atoms  and  a centrally  located  empty  spherical  region 
whose  e-xistence  is  essential  in  achieving  an  equilibrium  state 
for  the  cluster.  Some  modifications  were  made  to  calculate 
the  sphere  radii  of  the  overlapping  atomic  centers  self- 
consistently , and  also  to  include  a Madelung  correction  to 
atomic  potentials  due  to  the  ions  of  the  crystal  which  are  not 
represented  by  the  cluster  m.odel.  The  normalization  of  the 
partial-wave  components  and  acceleration  of  convergence  by 
a modified  Pratt  improvement  scheme  were  also  made  in  the  pro- 
graun.  From  this  model  bulk  properties  of  the  LiF  crystal  at 
equilibrium  were  obtained  in  very  good  agreement  with  experi- 
ments. To  conform  to  the  energy-band  structure  of  the  crystal, 
the  one-electron  energy  levels  of  the  cluster  were  linearly 
expanded.  The  F-center  proble.m  of  the  crystal  was  also  treated 
by  introducing  an  F-vacancy  at  an  F atomic  center  of  the 
OLi^F^  cluster,  and  the  consequent  change  of  the  electronic 
structures  and  the  optical  absorption  due  to  the  F-center  were 


studied. 
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The  results  show  that  the  SC?  Xj  a’j.tiple-scattering 
method  can  be  used  for  small  cluster  models  to  simulate  crystals 
(for  a mostly  covalent  crystal,  a GaP  crystal  is  treated  in 
a following  paper^)  and  to  study  point  defects.  Larger 
cluster  models  may  be  used  in  such  investigations  to  get  more 
accurate  quantitative  results. 

The  total  energies  of  clusters  made  of  heavy  atoms  are 
too  large  to  be  accurate  within  1 eV,  which  is  the  order  of 
magnitude  for  the  important  thermodynamic  parameters  such 
as  cohesive  energies  of  crystals  a.nd  formation  or  migration 
energies  of  point  defects.  To  reduce  the  magnitude  of  total 
energies  by  subtracting  common  terms,  we  can  freeze  core 
electrons  so  that  a nucleus  and  its  core  electrons  enter  as 
only  one  entity  in  calculating  various  terms  of  the  total 
energies.  Then  accurate  total  e.nergies  can  be  compared  to 
produce  the  thermodynamic  parameters  for  any  manageable-size 
clusters  of  any  atoms.  The  change  of  total  energy  by  changing 
the  atomic  positions  in  the  OLi^F^  cluster  was  not  accurate 
enough  to  incorporate  relaxation  effects  of  the  F-vacancy  in 
the  present  study;  this  situation  can  also  be  remedied  by 
adopting  the  proposed  scheme  of  calculating  the  total  energies. 
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TABLE  I 

Atomic  and  ionic  radii  of  Li  and  F,  in  a.u. 


Li 

F 

Atomic  radius  in  compounds 

2.740 

0.945 

Atomic  radius  in  tetrahedral 
covalent  bond^ 

— 

1.209 

Ionic  radius  of  Zachariasen 

1.285 

2.513 

* J.C.  Slater,  J.  Chem.  ?hys.  £1,  3199  (1964). 

^ L.  Pauling,  The  Nature  of  the  Chemical  Bond,  third  ed.,  Cornell 
Univ . Press,  New  York  (1960),  p.  246. 

From  C.  Kittle,  Introduction  to  Solid  State  Physics,  fourth  ed . , 
John  Wiley  & Sons,  New  York  (1971),  p.l29. 
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TABLE  II 


Results  of  the  SCF  Xa  cluster  model  calculation  for  the  LiF  molecule  with 
the  experimental  equilibrium  Internuclear  distance  (2.854  a.u.),  compared 


with  the  SCF  LCAO-MO  results  of  Brener  (ref.  29). 


Xa  Cluster 
Non-SCF 

Xa  Cluster 
SCF 

LCAO-MO 

SCF 

Radius  (a.u.) 


Li  center 

1.285 

1.351 

F center 

2.513 

2.643 

outer sphere 

3.326 

3.456 

Overlapping  scale 
factor 

1.331 

1.4 

One-electron  energies  (Ry) 

la  (F  Is  core) 

-49.168 

-48.207 

-52.252 

2a  (Li  Is  core) 

- 3.819 

- 3.658 

- 4.362 

3a 

- 2.498 

- 1.803 

- 2.639 

4a 

- 1.141 

- 0.477 

- 0.806 

llT 

- 1.151 

- 0.491 

- 0.697 

Constant  potential  in  the 
Interatomic  region  (Ry) 

- 0.440 

- 0.317 

Statistical  energies  (Ry) 

U,_,  (Nucleus-Nucleus) 

NN 

18.924 

18.924 

Ugjj(Electron-Nucleus) 

-571.021 

-554.683 

Ogg (Elec  tron-Elec  cron) 

140.617 

131.047 

Djj  (Exchange  correlation) 

- 25.550 

- 24.204 

(Kinetic) 

223.150 

214.694 

E.J  (Total  energy) 

213.879 

-214.223 

-212.763 

Virlal  coefficient  (-U  /E  ) 

T 

1.043 

1.002 

lonicity 

0.97 

0.92 
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TABLE  III 


Normalized  partial-wave  coefficients  (C^)  and  electronic  charge 
distributions  (Q)  of  the  valence  electron  states  of  the  LiF  molecular 
cluster  with  experimental  equilibrium  internuclear  distance. 


3a 

4a 

iTT 

C^:  F s 

0.993 

0.010 

L 

P 

-0.007 

0.974 

0.965 

Li  s 

0.061 

0.094 

— 

a P 

Out  s 

-0.059 

-0.038 

-0.135 

0.059 

0.048 

P 

-0.035 

0.115 

0.156 

d 

-0.004 

0.101 

0.055 

f 

-0.007 

0.003 

0.019 

Q^:  F 

0.987 

0.949 

0.931 

Li 

0.007 

0.027 

0.002 

® "Out"  stands  for  "extra-molecular". 

^ The  rest  of  the  electronic  charges  are  distributed  in  the  inter- 
atomic and  extramolecular  regions  of  the  cluster. 


TABLE  IV 


Results  of  the  SCF  St  cluster  model  calculations  for  the  LiF  molecule 


with  various  internuclear 

distances . 

Overlapping 

scale  factor 

is  1.4. 

Internuclear 

0.98  R * 

R 

1.02  R 

1.05  R 

distance 

e 

. e 

e 

e 

One-electron 

energies  (Ry) 

la  (F  Is  core) 

-48.215 

-48.207 

-48.202 

-48.195 

2o  (Li  Is  core) 

- 3.634 

- 3.653 

- 3.681 

- 3.713 

3a 

- 1.811 

- 1.803 

- 1.795 

- 1.784 

4a 

- 0.483 

- 0.477 

- 0.472 

- 0.464 

iTT 

- 0.500 

- 0.491 

t 

- 0.482 

- 0.470 

Constant  potential  in 

the  interatomic 

- 0.332 

- 0.317 

- 0.303 

- 0.284 

region  (Ry) 

Statistical  energies  (Ry) 

U,_.(Mucleus-Nucleus) 

19.3101 

18.9239 

18.5528 

18.0227 

U„(Eleccron-Nucleus) 

Cti^ 

-555.5654 

-554.6833  ■ 

-553.8552 

-552.6865 

U -(Electron-Electron) 

131.4323 

131.0467 

130.6818 

130.1655 

Oy  (Exchange  correladcn)-  24.2235 

- 24.2040 

- 24.1366 

- 24.1531 

(Kinetic) 

214.8266 

214.6943 

214.5844 

214.4404 

E^  (Total  energy) 

-214.2200 

-214.2225 

-214.2227 

■214.2208 

Virlal  coefficient 

1.0028 

1.0022 

1.0017 

1.0010 

(-U^/E^) 

lonlcity  (electronic 

0.9164 

0.9203 

0.9237 

0.9280 

charge) 

R^  indicates  the  experimental  equilibrium  internuclear  distance. 


"ifwiuvii'.'j 
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TABLE  V 


The  equilibriun  properties  of  the  LiF  molecule  as  calculated  by  the  SCF 
multiple-scattering  method,  compared  with  experimental  data. 


Calculation 


Experiment 


Internuclear  distance  (A) 


Vibrational  quantum 
energy,  (eV) 

Total  energy  (Ry) 

Cohesive  energy  (Ry) 


Dipole  moment  (Debye) 


1.53 

0.1197 

-214.223 

-0.539 

6.777 


1.51 

0.1124^ 

-215.004^ 

-0.440^ 

a 

6.6  + 0.3 


From  ref.  20. 

^ A.G.  Gaydon,  Dissociation  Energies,  rev.  ed.  (1953). 

^ 1 (Debye)  = 10  ^ (esu)  =*  1/4.303  (e-A)  ■ 0.5564  (a.u.) 
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TABLE  VI 

Results  of  the  SCF  Xot  aultlple-scattering  calculations  for  the  OLi^F^ 
cluster  with  various  nearest-neighbor  distances. 


Nearest-neighbor 

distance 

R • 
e 

1.05  R 

e 

1.074  R 

e 

1.1  R 1.15  R 

e e 

One-electron  energies  (Ry) 


Core 


F Is 

Li  Is 

-48.345 
- 3.624 

-48.333 
- 3.660 

-48.327 
- 3.677 

-48.318 
- 3.694 

-48.299 
- 3.725 

Valence  (Td  group) 
laj^  (F  2s) 

- 1.988 

- 1.959 

- 1.946 

- 1.933 

- 1.908 

lt2  (F  2s) 

- 1.973 

- 1.948 

. - 1.937 

- 1.925 

- 1.903 

2a^  (r  2?) 

- 0.716 

- 0.682 

- 0.668 

- 0.652 

- 0.622 

le  (F  2p) 

- 0.675 

- 0.646 

- 0.634 

- 0.621 

- 0.595 

2t2  (F  2p) 

- 0.672 

- 0.646 

- 0.635 

- 0.622 

- 0.597 

Itj^  (F  2p) 

- 0.632 

- 0.612 

- 0.603 

- 0.593 

- 0.573 

3t2  (F  2?) 

- 0.620 

- 0.602 

- 0.594 

- 0.584 

- 0.566 

Constant  potential 
in  the  inter- 

- 0.363 

- 0.329 

- 0.315 

- 0.300 

- 0.274 

atomic  region 
Statistical  Energies (Ry) 


U„(Madelung  correction) 

- 0.3857 

- 0.3905 

- 0.3916 

- 0.3921 

- 0.3914 

O.r.,  (Nuc  leus-Nuc  leus  ) 

404.4785 

385.2187 

376.6104 

367.7037 

351.7214 

U_„ (Electron-Nucleus) 

C*Cl 

-2873.1807 

-2834.3291 

-2817.0615 

-2799.1934 

-2767.1826 

Ugg (Electron-Elec  cron) 

850.9722 

831.9631 

823.5098 

814.7700 

799.1123 

(Exchange  correlation) 

- 96.4353 

- 96.3403 

- 96.3048 

- 96.2679 

- 96.2084 

Uj^  (Kinetic) 

857.3088 

856.5955 

856.3496 

865.0928 

855.6870 

Ey  (Total  energy) 

-857.2419 

-857.2825 

-857.2331 

-857.2817 

-857.2612 

Virial  coefficient 

1.0001 

0.9992 

0.9989 

0.9986 

0.9982 

lonicity  (electronic 

0.8812 

0.8935 

0.8988 

0.9041 

0.9133 

charge) 


is 


the  experimental  nearest-neighbor  distance  at  equilibrium. 
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TABLE  VII 


The  equilibrium  properties  of  the  LiF  crystal  as  calculated 
by  the  SC?  Xa  multiple-scattering  method  for  the  OLi^F^ 
cluster  model,  compared  with  experimental  data. 


OLi4F4  cluster 
calculation 

LiF  crystal 
experiment 

5 

Radii  of  atomic  centers  (A) 

Li 

1.086  (0.952)^ 

F 

1.376  (1.207) 

Nearest-neichbor 

2.159 

^ b 

distance  (.^i) 

2.01 

Bulk  modulus  (dyne/cm^ ) 

6.66  X 10' * 

6.98  X 10"*^ 

Energies  per  molecule  (Ry) 

■ 

Total  energy 

-214.322 

Cohesive  energy 

0.638 

0.644“^ 

lonicity  (electronic  charge) 

0.999 

The  numbers  outside  the  parentheses  are  the  radii  of  the 
atomic  spheres  determ, ined  through  SCF  calculations  by  satis- 
fying the  virial  theorem.  The  overlapping  scale  factor  is 
1.14.  The  numbers  inside  the  parentheses  are  the  radii  of 
the  touching  atomic  spheres  which  are  reduced  from  the 
overlapping  atomic  spheres  so  that  the  overlapping  scale 
factor  is  u.nity. 

From  ref.  23. 

See  text. 


j 

i 
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TABLE  VIII 


Normalized  partial-wave  coefficients-  (C^)  and  electronic  charge  dis- 
tributions  (Q)  of  the  St^,  3a^,  and  states  of  the  OLi^F^  cluster. 


^^2 

"^2 

F s 

-0.026 

0.278 

-0.139 

P 

0.862,  0.465 

-0.013 

-0.126,  0.033 

Li  s 

-0.021 

-0.237 

-0.120 

P 

-0.044,  0.044 

0.035 

0.118,  0.004 

0 s 

— 

0.064 

— 

P 

-0.043 

— 

0.009 

Out^  s 

— 

0.766 

— 

P 

0.085 

— 

0.890 

d 

-0.047 

— 

0.050 

f 

0.027  

0.046 

-0.022 

Q 

F 

0.959 

0.077 

0.036 

Li 

0.004 

0.057 

0.023 

0 

0.002 

0.004 

0.0 

In'* 

0.024 

0.272 

0.141 

^ a 

Out 

0.010 

0.589 

0.794 

a,, 


Out"  stands  for  "extramolecular". 
"In"  represents  the  interatomic 


rsr> 
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O 

TABLE  IX 


Compatibility  relations  between  the  representations  of  the  Td 
and  point  groups. 
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TABLE  X 


Results  of  the  SCF  Xa  multiple-scattering  calculations  for  the  VacF 
(F-vacancy)  cluster  in  various  configurations. 


Cluster 

(VacF)'*’^ 

(VacF) 

(VacF) 

spin-polarized 

configuration 

spin- 

restricted 

spin 

t spin 

One-electron 
energies  (Ry) 

Core  states 


F Is 

-48.369 

-48.341 

-48.340  -48. 

Li  Is 

- 3.811 

- 3.713 

- 

3.707  - 3. 

ALi  Is 

- 3.672 

- 3.672 

- 

3.674  - 3. 

Valence  states 

la-  (F  2s) 

- 1.982 

- 1.962 

- 

1.956  - 1, 

le^  (F  2s) 

- 1.975 

- 1.955 

- 

1.949  - 1, 

2ai  (F  2p) 

- 0.697 

- 0.674 

- 

0.671  - 0. 

2e  (F  2?) 

- 0.672 

- 0.652 

- 

0.648  - 0. 

3aj^  (F  2?) 

- 0.663' 

- 0,644 

- 

0.640  - 0. 

3e  (F  2p) 

-• 0.655 

- 0.635 

- 

0.631  - 0. 

la_  (F  2p) 

- 0.641 

- 0.620 

- 

0.616  - 0. 

4e  (F  2?) 

- 0.632 

- 0.612 

- 

0.608  - 0. 

Excited  states 

4ai 

- 0.188 

- 0.133 

- 

0.225  - 0. 

5ai 

- 0.031 

- 0.053 

- 

0.066  - 0. 

5e 

- 0.013 

- 0.034 

- 

0.047  - 0.1 

Constant  potential  in 

the  interatomic 

- 0.327 

- 0.337 

- 

0,354  - 0. 

region 

Statistical  energies  (Ry) 

U„  (Madelung  correction) 

M 

- 0.3555 

- 0.3760 

- 0.3797 

Ujjj^(Nucleus-Nucleus) 

245.0178 

245.0173 

245.0178 

0_..  (Electron-Nucleus) 

EN 

-2057.0527 

-2071.1079 

-2071.4536 

Ugj.  (Electron-Electron) 

573.6360 

587.4375 

587.7170 

(Exchange  correlation) 

- 75.5144 

- 75.6286 

- 75.6914 

(Kinetic) 

656.6375 

656.8943 

655.9946 

(Total  energy) 

-657.5808 

-657.7620 

-657.7947 

lonicty  (electronic  charge) 

VacF 

- 0.022 

- 0.769 

- 0.794 

F 

- 0.884 

- 0.914 

- 0.912 

Li 

0.924 

0.871 

0.873 

ALI 

0.903 

0.898 

0.398 
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TABLE  Xlb 

Nonnallzed  parial-wave  coefficients  (C^)  and  electronic  charge  distrib- 

la 

utions  (Q)  of  the  ^a^^  and  states  for  both  spins  in  the  equilibriun 
VacF  cluster. 


F s 

0.125 

0.207 

0.223 

-0.149 

Pi 

0.020 

0.012 

-0.055 

0.076 

P2 

-0.158 

-0.091 

0.081 

-0.101 

Li  s 

-0.216 

-0.243 

-0.096 

-0.002 

Pi 

-0.191 

-0.093 

0.119 

-0.116 

P2 

0.040 

0.045 

-0.035 

0.070 

ALi  s 

-0.027 

-0.076 

-0.150 

0.134 

P 

0.011 

0.004 

-0.051 

0.058 

VacF  s 

0.660 

0.390 

-0.153 

0.206 

P 

0.0 

-0.007 

-0.036 

0.032 

0 s 

0.110 

0.100 

0.026 

-0.004 

P 

0.043 

0.022 

-0.016 

0.018 

Out^  s 

0.306 

0.602 

0.636 

-0.407 

P 

0.242 

0.215 

-0.458 

0.720 

d 

-0.138 

-0.077 

0.087 

-0.062 

F (3) 

0.041 

0.051 

0.059 

0.033 

Li  (3) 

0.085 

0.074 

0.025 

0.018 

ALi 

0.001 

0.006 

0.025 

0.021 

VacF 

0.436 

0.152 

0.025 

0.044 

0 

0.014 

0.010 

0.001 

0.0 

Inb 

0.252 

0.294 

0.243 

0.191 

Out® 

0.171 

0.141 

0.622 

0-.637 

"Out"  stands  for  "extramolecular" 
^ "In" 


represents  the  interatomic  region  of  the  cluster. 


FIGURE  CAPTIONS 


Fig.  1 
Fig.  2a 


Fig.  2b 
Fig.  2c 
Fig.  3 


Cluster  model  for  a LiF  molecule. 

Contour  map  of  the  3a  wave  function  of  the  LiF  molecule. 
The  spheres  of  the  atomic  centers  and  the  outersphere 
are  also  shown.  The  contour  lables  indicate  the 
wave  function  values  whose  m.agnitudes  are  equally 
spaced  in  logarith.mic  scale  between  extrema.  The 
coordinates  are  in  a-u.,  centered  at  the  center  of 
the  outersphere. 

Controu  map  of  the  4a  wave  function  of  the  LiF  molecule. 
Controu  map  of  the  It  wave  function  of  the  LiF  molecule. 
Calculated  total  energy  of  the  LiF  molecular  cluster 
vs.  internuclear  distance  in  units  of  the  experimental 
equilibrium  internuclear  distance.  The  vertical  bars 
indicate  the  ranges  of  the  total  energies  at  the 
later  stages  of  iterations.  The  dotted  curve  extends 
a parabolic  fit. 


Fig.  4 The  OLi^F^  cluster  model  for  a LiF  crystal. 

Fig.  5 Calculated  total  energy  of  the  OLi^F^^  cluster  vs.  lattice 
constant  in  units  of  the  experi.mental  equilibrium 
lattice  constant  of  the  LiF  crystal.  The  vertical 
bars  indicate  the  ranges  of  the  total  energies  at  the 
later  stages  of  iterations. 
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Fig . 6 


Fig.  7a 


Fig.  7b 


Fig.  7c 

Fig.  8 


One-electron  energy  levels  of  the  OLi^F^  cluster 
in  the  ground  state  configuration  (a) , and  in  the 

transition  state  configuration  (b) . represents 

the  constant  potential  in  the  interatomic  region 
of  the  cluster.  The  numbers  in  parentheses  in- 
dicates the  occupation  numbers  of  the  band  edge 
states;  all  the  levels  below  them  are  completely 
filled,  and  those  above  them  are  empty. 

Contour  map  of  the  3t2  wave  function  of  the  OLi^F^ 
cluster.  It  is  plotted  in  the  llo  plane  (refer 
to  Fig.  4)  which  bisects  two  Li  and  two  F atoms, 
and  the  broken-line  circles  are  the  projections 
of  the  atomic  centers  (including  the  0 center)  and 
the  outersphere  of  the  cluster.  The  contour 
labels  indicate  the  wave  function  values  whose 
magnitudes  are  equally  spaced  in  logarithmic  scale 
between  extrema. 

Contour  map  of  the  33^^  wave  function  of  the  OLi^F^ 
cluster. 

Contour  mao  of  the  4t.,  wave  function  of  the  OLi^F. 
cluster . 

Linearly  expanded  one-electron  energy  levels  of  the 
OLi^F^  cluster  in  (a)  to  fit  the  bandgap  of  the  LiF 
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crystal  whose  energy  bands  are  shown  in  (b)  as 
calculated  by  Brener  (ref.  29),. 

Fig.  9 One-electron  energy  levels  of  the  VacF  (F-vacancy) 
cluster  in  various  configurations  (b ,c , and  d) , com- 
pared with  those  of  the  OLi^F^  cluster  (a) . 

represents  the  constant  potential  in  the  interatomic 
region  of  the  cluster.  The  names  of  the  valence 
electron  states  for  the  VacF  cluster  in  all  the  con- 
figurations are  the  same  in  the  same  order.  Shaded 
levels  are  the  highest  occupied  ones  with  their 
occupation  numbers  being  indicated  by  the  portions 
of  the  shades. 

Fig.  lOaj^  Contour  map  of  the  4e^  wave  function  of  the  VacF 

cluster.  It  is  plotted  in  the  plane  which  bisects 

the  F-vacancy,  .A.Li  atom,  and  each  one  of  the  Li  and 

F atoms.  (It  is  the  same  plane  as  the  one  used  for 

the  contour  maos  of  the  wave  functions  of  the  OLi.F.  i 

* 4 4 

electron  states  — Figs.  7a,  lb,  and7c  — while  the 
upper  right  F atomic  certter  is  now  replaced  by  the 
F-vacancy.)  The  broken-line  circles  are  the  pro- 
jections of  the  atomic  centers  and  the  outersphere 
of  the  cluster.  The  contour  labels  indcate  the 
wave  function  values  whose  magnitudes  are  equally 
spaced  in  logarithmic  scale  between  extrema. 
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Fig. 

10a2 

Contour 

map 

of 

the 

• -V 

4e  wave 

function  of 

the  ' 

VacF  1 

cluster . 

Fig. 

lObj^ 

Contour 

map 

of 

the 

wave 

function 

of 

the 

VacF 

cluster . 

Fig . 

10b2 

Contour 

map 

of 

the 

wave 

function 

of 

the 

VacF 

cluster . 

Fig. 

lOCj^ 

Contour 

map 

of 

the 

5^1 

wave 

function 

of 

the 

VacF 

cluster . 

Fig. 

IOC2 

Contour 

map 

of 

the 

Sa^^ 

wave 

function 

of 

the 

VacF 

cluster . 

Fig. 

lOd^ 

Contour 

map 

of 

the 

5e^ 

wave 

function 

of 

the 

VacF 

cluster . 

Fig. 

10d2 

Contour 

map 

of 

the 

Se'*' 

wave 

function 

of 

the 

VacF 

cluster . 

Fig . 

11. 

Potentials  : 

for 

spin 

i-up 

and 

spin-down 

electn 

cns  in  . ... 

the  F-vacancy  region  of  the  VacF  cluster.  The 

horizontal  lines  indicate  the  constant  potentials 
in  the  interatomic  region  of  the  cluster. 
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ABSTRACT 


A PbTe  crystal  is  modeled  by  an  OPb^Te^  cluster,  comprised  of 
four  tetrahedrally  placed  atoms  for  each  of  Li  and  F,  and  a central 
empty  spherical  region  called  0.  Then  vacancies,  interstitials,  and 
Frenkel  defects  of  either  Pb  or  Te  are  introduced  in  the  cluster, 
resulting  in  defect  clusters  which  simulate  the  PbTe  crystals  with 
such  point  defects.  The  various  clusters  are  investigated  by  the  SCF 
Xa  multiple-scattering  method  of  Johnson,  with  the  mass-velocity  and 
Darwin  relativistic  corrections.  The  results  that  ?b  and  Te  vacancies 
produce  two  free  holes  and  two  free  electrons  agree  with  Parada  and 
Pratt's  results  of  energy-band  calculations  through  the  Green's 
function  method.  .\  Pb  interstitial  and  a Pb  Frenkel  defect  are  shown 
to  give  four  and  two  free  electrons  respectively.  A Te  interstitial 
seems  to  generate  two  free  holes  while  a Te  Frenkel  defect  does  not  seem 
to  contribute  any  to  the  conductivity.  Anti-bonding  mechanisms  are 
found  to  play  very  important  roles  in  the  interstitial-containing  com- 
plexes by  shifting  some  valence  electron  levels  of  a pus  crystal  up  into 


the  conduction  band. 
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1.  INTRODUCTION 


PbTe,  PbSe,  and  PbS  fora  a group  a lead  chalcogenldes.  Group  IV-VI 


semiconductors  of  NaCl  structure.  They  have  been  the  subject  of  extensive 
investigations  because  of  their  technologically  Important  optical  and 
thermoelectric  properties.  These  properties  are  closely  related  to  the 
free  carriers  produced  by  the  intrinsic  point  defects  which  are  controlled 
by  the  use  of  high  ranges  of  non-stoichiometry. 

Experiments  on  the  effect  of  an  intrinsic  point  defect  on  the  con- 
ductivity of  lead  chalcogenides  lead  to  ambiguous  conclusions  since 
there  usually  are  sef/eral  kinds  of  intrinsic  point  defects  present  simul- 
taneously in  one  sample.  The  are  reviewed  by  Abrikosov,  et  al.^.  Ravlch, 

2 3 

et  al.  and  Davin.  In  this  paper  it  is  attempted  to  solve  the  problem 
theoretically  for  vacancies,  interstitials,  and  Frenkel  defects  of  both 
Pb  and  Te  in  PbTe. 

Either  p-type  or  n-type  semiconductor  can  be  made  of  PbTe  by  adding 
excess  tellurium  or  lead.  In  both  cases  the  free  carriers  are  char- 
acterized by  high  mobilities,  and  they  can  not  be  frozen  out  even  at 

4 5 

very  low  temperatures.  ’ With  such  experimental  evidence,  it  has  been 
assumed  that  there  are  no  bound  states  in  the  bandgap.  Very  recently, 
Bauer,  et  al.^  showed  conclusively  that  their  measurements  of  the  Hall 
effect  and  the  magnetoresistance  rule  out  the  existance  of  deep  impurity 
levels  as  well  as  of  levels  originating  from  intrinsic  point  defects 
within  the  bandgap  in  PbTe.  If  there  were  any  bound  states,  in  Che  band- 
gap,  it  would  have  been  very  difficult  to  exhibit  them  in  calculations 
because  the  bandgap  is  extremely  narrow.  It  is  only  0.19  eV  at  4.2*K, 


Increasing  to  0.32  eV  at  300*K.^ 

The  problens  of  Pb  and  Te  vacancies  in  PbTe  were  first  solved  by 
8“10 

Parada  and  Pratt  by  the  Green's  function  method  of  Koster  and 

Slater^^  in  which  a point  defect  is  treated  by  a localized  pertiirbing 

potential  and  the  pertxirbed  wave  functions  of  electrons  are  expanded 

in  terms  of  the  Wannier  functions  of  the  unperturbed  lattice.  The 

calculations  are  quite  complicated  because  the  Bloch  functions  as  well 

as  the  eigen-energies  have  to  be  known  for  all  k vectors  of  the  entire 

Brillouin  zone  in  order  to  compute  the  Wannier  functions.  Therefore, 

the  calculations  were  carried  out  without  self-consistency  using  the 

potentials  of  natural  atoms. 

• 12 

In  recent  years,  Johnson  has  developed  the  SCP  Xa  multiple-scattering 

method  for  molecular  clusters.  The  method  is  very  efficient  and 

accurate  in  comparison  with  the  traditional  methods  such  as  the  LCAO 

method;  it  has  been  successfully  applied  to  many  polyatomic  molecules. 

It  is  also  possible  to  treat  crystals  by  the  method  once  suitable  cluster 

models  are  chosen  for  the  crystals.  Vacancies  in  PbTe  were  considered 
13 

by  Hamstreet  using  a 27  atom  cluster  for  the  crystal.  In  this  paper, 

the  smallest  possible  neutral  cluster  of  eight  atoms  and  a central  empty 

spherical  region  is  employed  in  simulating  the  PbTe  crystal,  and  all  the 

intrinsic  point  defects  are  studied  within  the  same  cluster  format.  This 

cluster  model  is  the  one  which  was  adopted  for  a LiF  crystal  in  a previous 

14 

paper  by  the  present  authors. 

In  Section  II,  the  cluster  model  calculation  for  a perfect  PbTe 
crystal  is  described;  and  subsequently  in  Sections  III-V,  vacancies, 
interstitials,  and  Frenkel  defects  of  Pb  and  Te  are  treated. 
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II.  OPb.Te,  CLUSTER  MODEL  CALCULATION 
4 4 

The  OPb^Te^  cluster  model  for  a PbTe  crystal  is  schematically 

shown  in  Fig.  1,  corresponding  to  a cubic  section  of  a PbTe  crystal. 

The  symmetry  of  the  cluster  now  belongs  to  the  T,  point  group  reduced 

a 

from  the  Oj^  for  the  crystal.  The  central  empty  spherical  region,  0, 

is  inserted  so  that  the  scattered  waves  from  the  atomic  centers  are  more 

accurately  described  in  the  central  region  of  the  cluster  by  the  use  of 

a spherically  symmetric,  r-dependent,  potential  rather  than  a volume- 

averaged  one.  The  Importance  of  the  0 center,  was  illustrated  in  the 
14 

previous  paper  for  the  LiF  crystal  in  achieving  an  equilibrium  state 

for  the  OLi^F^  cluster.  The  nearest-neighbor  distance  of  the  cluster 

. • 

model  was  taken  from  the  experimental  equilibrium  value  3.226  A of  the  crystal. 

In  the  OPb^Te^  cluster,  Pb  and  Te  atoms  are  equally  treated  and 
every  atom  is  equivalent  to  the  other  atoms  of  the  same  species  as  they 
should  be  in  the  PbTe  crystal.  However,  each  atom  has  three  nearest- 
neighbors  in  the  cluster  while  it  has  six  in  the  crystal.  Therefore, 
there  will  be  some  rearrangement  in  the  covalent-bonding  characteristics, 
but  the  cluster  model  is  quite  adequate  in  describing  the  ionic-bonding 
characteristics  with  the  use  of  Madelung  correction  terms  to  the  atomic 
potentials  from  the  ions  in  the  remainder  of  the  crystal  which  are  not 
represented  in  the  cluster  model  (refer  to  Ref.  14).  Due  to  the  neutrality 
of  the  cluster  all  the  electrons  can  be  contained  without  the  need  of  [ 

fi 

using  a Watson  sphere. 

The  atomic  radii  of  Pb  and  Te  in  compounds  are  1.80  and  1.40  A 
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respectively.^^  However,  the  atooic  potentials  calculated  self- 
consistently  for  the  Pb  and  Te  atomic  centers  in  the  cluster  model  are 
almost  the  same  at  the  midpoint  between  the  two  adjacent  atoms.  There- 
fore, both  the  atomic  radii  were  set  equal  in  the  subsequent  iterations. 

14 

In  the  case  for  the  OLi^F^  cluster,  the  overlapping-scale  factor 
of  the  atomic  centers  was  determined  such  that  the  virial  coefficient^^ 
becomes  unity.  The  virial  coefficient  calculated  for  the  OPb^Te^  cluster 
with  touching  atomic  centers  is  already  very  close  to  unity  with  a de- 
viation of  less  than  0.01*  due  to  the  enormously  deep  core  energy  levels 
(for  example,  the  eigen-energy  of  Pb  Is  core  electron  state  is  about 
-6000  Ry).  Because  the  virial  coefficient  hardly  varies  by  changing  the 
degree  of  the  overlap  for  the  atomic  centers,  the  overlapping  scale 

factor  of  the  OP.Te,  cluster  is  set  at  1.14  as  we  found  for  the  OLi.F. 

4 4 4 4 

cluster. 

The  exchange  parameter  oh  for  Pb  and  Te  were  djtained  through  inter- 

18 

polation  from  the  Schwarz's  atomic  a's.  For  the  interatomic  and  extra- 

molecular  regions,  a weighted  average  cc  was  computed  with  respect  to 

the  number  of  valence  electrons  of  the  atoms.  In  the  OPb.Te,  cluster 

all  the  a values  are  approximately  0.7  which  is  close  to  the  lower  limit. 

It  was  found  that  relativistic  effects  are  important  in  getting 

19  20 

the  correct  energy-band  structure  of  the  PbTe  crystal.  ' Therefore, 
we  Included  the  mass-velocity  and  Darwin  relativistic  corrections  to  the 
SCF  Xa  multiple-scattering  method  for  the  OPb^Te^  cluster,  (the  formalism 
is  Illustrated  in  Appendix) . The  spin-orbit  coupling  term  was  not  in- 


6 

corporaced  because  of  the  synmetry  restrictions.  For  a given  potential, 
the  relativistic  corrections  tend  to  lower  the  electron  energy  levels; 
however,  the  relativistic  and  non-relativistic  energy-level  structures 
become  Increasingly  similar  as  self-consistent  iterations  proceed. 

The  calculated  electron  states  of  the  OPb.Te,  cluster  are 

listed  in  Table  I with  their  eigen-energies,  normalized  partial-wave 

coefficients,  and  electronic  charge  distributions.  The  electron  energy 

levels  of  the  cluster  are  shown  in  Fig.  2(a),  compared  with  the  energy- 

band  structure  of  the  PbTe  crystal  in  Fig.  2(b)  which  was  obtained  by 

the  pseudopotential  method  using  five  parameters  to  fit  experimental 

21 

data  by  Lin  and  Kleinman. 

The  two  lowest  levels  of  the  cluster,  la^^  and  lt2,  correspond  to 
the  lowest  crystal  energy-band  of  mostly  Te  s character;  and  the  next 
two  levels,  2t2  and  Za^^,  correspond  to  the  second  lowest  crystal  energy- 
band  of  mostly  Pb  s character.  These  two  energy-bands  will  be  referred 
to  as  Te 3 and  Pb  s bands  for  convenience. 

The  next  five  levels,  le  through  ^^2,  form  the  top  valence  band 
which  is  primarily  made  of  Te  p basis  states  but  they  also  contain 
significant  components  of  Pb  p basis  states.  The  rest  of  the  levels  from 
5C2  and  up  constitute  the  conduction  band  of  the  crystal  as  they  are 
made  of  mostly  Pb  p basis  states  and  some  Te  p basis  states.  The 
conduction  band  states  have  large  amounts  of  electronic  charge  in  the. 
interatomic  region. 

As  we  see  from  Fig.  2,  each  energy-band  of  the  crystal  matches  well  to 
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a sec  of  cluster  energy  levels  of  similar  characterlsclcs . The  bandgap 
of  the  cluster  as  the  difference  in  eigen-energies  between  the 
and  5C2  levels,  is  0.19  Ry  which  is  about  ten  times  as  large  as  the 
bandgap  of  the  crystal  at  L point  of  the  Brillouin  zone.  The  band- 
gap  of  the  cluster,  rather,  corresponds  to  an  average  of  the  energy  gap 
between  the  conduction  and  the  valence  bands  throughout  the  Brillouin 
zone. 

The  non-relativistic  total  energy  of  the  cluster  is  -209,088,7  Ry 
self -consistently , while  the  sum  of  the  RHF  total  energies  for  free 
atoms  is  -209,086.1  Ry.  This  gives  2.6  Ry  for  the  cohesive  energy  of 
the  cluster;  however,  it  should  be  considered  to  give  only  an  indication 
of  cohesion  since  the  total  energies  are  too  large  to  be  accurate  on  the 
order  of  0.01  Ry.  The  inclusion  of  the  relativistic  correction 

terms  drastically  lowers  the  total  energy  of  the  cluster  by  4671.7  Ry 
from  the  non-relativistic  value  while  the  corresponding  relativistic 
RHF  total  energies  of  free  atoms  are  not  available.  It  is  impossible 
to  determine  the  equilibrium  state  of  the  cluster  or  formation  energies 
of  defects  by  comparing  the  total  energies  of  clusters  which,  for  such 
purposes,  must  be  calculated  accurately  to  the  order  of  0.1  Ry. 

The  ionicty  of  the  OPb^Te^  cluster  indicates  that  about  0.75 
electron  is  transferred  from  a ?b  to  a Te  atom.  It  was  calculated  by 
assigning  the  electronic  charges  in  the  interatomic,  extramolecular , 
and  the  0 center  regions  to  the  atomic  centers  in  proportion  to  the 
valence  electronic  charges  which  are  already  in  the  atomic  centers 
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III.  Pb  AND  Te  VACANCIES 

A Pb  or  Te  vacancy  cluster  was  made  by  vacating  a Pb  or  Te  neutral 

atom  from  the  OPb.Te,  cluster.  This  results  in  reduction  of  the  cluster 
4 4 

symmetry  from  the  T,  to  C,  point  group.  (In  the  crystal  the  symmetry 

a jv 

is  reduced  from  the  0,  to  T.  point  group.)  The  exchange  a parameters 

u G 

are  the  same  as  the  OPb^Te^  cluster  except  in  the  vacancy  centers  where 
Che  average  a is  used.  The  same  Madelung  correction  terms  were  kept 
but  now  xjith  the  final  ionicity  found  with  the  OPb^Te^  cluster. 

The  energy  levels  of  the  Pb  and  Te  vacancy  clusters  are  shown  in 
Figs.  3(a)  and  3(c)  respectively,  compared  with  the  levels  of  the 
OPb^Te^  cluster  in  Fig.  3(b). 

We  will  discuss  the  results  6f  the  Pb  vacancy  cluster  first.  The 
lowest  three  levels  of  Fig.  3(a),  1®^^,  2a^,  and  le,  constitute  a perturbed 
Te  s band  due  to  the  Pb  vacancy.  Likewise,  the  next  two  levels,  2e  and 
3aj^,  constitute  a perturbed  Pb  s band.  The  set  of  the  levels  from  3e 
through  6e  make  up  the  top  valence  band,  and  the  levels  of  7a^  and  up 
make  up  the  conduction  band  of  the  perturbed  cluster. 

The  topmost  valence  electron  level  is  the  four-fold  6e  state,  and 
it  is  occupied  by  only  two  electrons;  therefore,  the  valence  band  has 
two  less  electrons  than  required  if  it  were  to  be  completely  filled. 

We  notice  that  an  a^  level  is  missing  from  the  Pb  s band  as  a Pb  is 
removed  from  the  OPb^Te^  cluster  by  considering  the  compatibility  re- 

lationships  for  the  representations  between  the  T^  and  point  groups.  J 

In  terms  of  a perturbative  treatment,  a doubly  degenerate  atomic  Pb  s ] 
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level  is  so  much  perturbed  by  the  vacancy  potential  that  Is  is  driven 
up  into  the  conduction  band.  Now  the  Pb  vacancy  cluster  has  two  less 
electron  states,  counting  degeneracy,  in  the  valence  bands  and  four 
less  valence  electrons  due  to  the  missing  Pb  atom  than  the  OPb^Te^ 
cluster;  therefore,  it  results  in  two  free  holes  in  the  valence  band. 

The  extremely  narrow  bandgap  of  PbTe  is  located  very  close  to  the 
constant  potential  in  the  interatomic  region  of  the  OPb^Te^  cluster  as 
we  see  in  Fig.  2.  Figure  3 now  shows  that  the  constant  potential  hardly 
varies  when  the  Pb  vacancy  (or  Te  vacancy  for  that  matter)  is  introduced 
in  the  OPb^Te^  cluster,  and  it  still  Indicates  approximately  the  position 
of  the  bandgap  of  the  crystal.  The  large  gap  between  the  6e  and  7a^ 
levels  further  assures  that  the  crystal  bandgap  must  be  between  them. 

From  the  partial-wave  analysis  it  is  seen  that  the  6e  state  is  indeed 
a valence-band  state  whose  characteristics  are  very  similar  to  those 
of  the  top  valence  electron  state,  4t2,  of  the  OPb^Te^  cluster. 

Now  we  will  discuss  the  case  for  the  Te  vacancy.  The  two  lowest  levels, 
laj^  and  le,  constitute  a perturbed  Te  s band.  This  band  is  lacking  one 
state  compared  with  the  Te  s band  of  the  OPb^Te^  cluster  accounting 
for  the  5s  level  of  the  missing  Te  atom.  The  next  three  levels,  2a^, 

2e,  and  3aj^,  constitute  a perturbed  Pb  s band. 

The  set  of  next  six  levexs,  from  3e  to  Sa^^,  makes  up  a perturbed 
top  valence  band  which  now  accommodates  only  18  electrons  while  there 
are  24  electron  states  in  the  corresponding  valence  band  of  the  OPb^Te^ 
clsuter.  The  effect  of  the  perturbation  by  the  Te  vacancy  is  to  make 
an  a^  and  e levels  missing  from  the  top  valence  band.  Skipping  the 
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6aj^  and  6e  levels  for  the  moment,  the  top  two  levels,  7e  and  7a^,  belong 
to  the  conduction  band  as  their  wave  functions  have  similar  characteristics 
as  the  wave  functions  of  the  conduction  band  levels  of  the  OPb.Te,  cluster. 

In  order  to  illustrate  the  fa^^  and  6e  states,  we  list  their  partial- 
wave  coefficients  and  electronic  charge  distributions  in  Table  II.  Both 
the  states  have  large  partial-wave  coefficients  for  the  Pb  p basis 
functions  and  considerable  amount  of  electronic  charges  in  the  inter- 
atomic region;  such  properties  are  characteristic  of  the  conduction  band 
states  of  PbTe.  They  also  lie  above  the  constant  potential  which  is  in- 
dicative of  the  location  of  the  crystal  bandgap.  Furthermore,  it  is 
known  experimentally  that  there  are  no  bound  electron  states  in  the 
bandgap  of  the  crystal;  therefore,  the  two  states  belong  to  the  con- 
duction band  of  the  PbTe  crystal. 

The  6aj^  and  6e  states  also  have  significant  contributions  from  Te 
vacancy  basis  functions.  .The  state  has  particularly  large  coefficient 
for  a Te  vacancy  s basis  function  so  that  it  can  be  said  to  be  the 
electron  state  of  the  greatly  perturbed  Te  5s  atomic  orbital  which  was 
stated  earlier  as  missing  from  the  perturbed  Te  s band.  On  the  other 
hand,  the  6e  state  is  the  perturbed  electron  state  which  has  been  driven 
out  of  the  top  valence  band  of  primarily  Te  p character.  The  ocher 
perturbed  a^^  level,  which  has  also  been  driven  out  of  the  top  valence 
band,  would  be  located  further  above  the  7a^  level. 

The  6aj^  conduction  band  level  is  the  highest  filled  level  occupied 
by  two  electrons.  The  situation  can  be  expalined  by  the  following.  Due 
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to  a missing  Te  atom  there  are  six  fewer  electrons  and  eight  fewer 
electron  levels,  counting  degeneracy,  in  the  valence  bands  of  the  Te 
vacancy  cluster.  Consequently  we  have  two  extra  electrons  left  over 
after  all  the  valence  band  levels  are  filled.  These  two  electrons  go  to 
the  lowest  conduction  band  level,  Sa^^,  resulting  in  two  free  electrons 
in  the  PbTe  crystal. 

Our  conclusions  on  the  effects  of  the  Pb  and  Te  vacancies  on  the 

conductivity  of  PbTe  are  in  agreement  with  the  results  of  perturbative 

8—10 

energy-band  calculations  of  Parada  and  Pratt  by  the  Green's  function 

method.  In  their  calculations  for  the  Te  vacancy,  the  perturbed  Te  s 

level  is  slightly  higher  than  the  perturbed  Te  p levels  while  the  present 

cluster-model  calculations  show  otherwise.  In  any  case  such  change  of 

ordering  does  not  present  and  change  in  physical  meaning  because  they 

are  all  in  the  conduction  band.  In  the  Green's  function  method,  all 

the  Bloch  states  of  the  crystal  are  considered  unperturbed  except 

producing  some  localized  defect  electron  states,  while  in  the  cluster 

method  all  the  states  are  perturbed  by  the  defect.  The  multiple-scattering 

13 

calculations  for  27-atom  clusters  by  Hamstreet  also  give  the  same 
results  regarding  to  the  effects  of  the  vacancies  on  the  conductivites 
of  the  PbTe  crystal. 
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IV.  Pb  AND  Te  INTERSTITIALS 

22 

In  PbTe  the  Pb  interstitial  is  still  important,  although  the 

primary  intrinsic  defect  species  is  the  Te  vacancy  when  Che  material 

23 

contains  excess  lead.  It  seems  that  the  Pb  interstitials  can  be 
rather  easily  accommodated  by  the  crystal  due  to  their  small  ionic 
size,  but  the  Te  ions  are  too  large  to  be  interstitials.  However,  for 
competeness,  both  the  Pb  and  Te  interstitials  will  be  considered  here. 

The  interstitial  clusters  are  formed  very  conveniently  by  putting 
either  Pb  or  Te  atom  at  the  place  of  the  central  0 empty  center  of  the 
OPb^Te^  cluster.  The  resulting  interstitial  clusters  retain  the  T^ 
group  symmetry  which  is  correct  for  interstitials  in  the  crystal.  The 
electron  energy  levels  of  the  ?b  and  Te  interstitial  clusters  are  shown 
in  Fig.  4,  compared  with  chose  of  the  OPb^Te^  cluster. 

Let  us  consider  the  Pb  interstitial  case  first.  From  the  partial- 
wave  analysis  it  can  be  shown  that  all  the  e,  t^,  and  t^  states  of  the 
Pb  interstitial  cluster  are  perturbed  electron  states  from  the  ones  with 
the  same  denominations  in  the  OPb^Te^  cluster.  Their  eigen-energies  have 
been  lowered  by  almost  Che  same  amount  with  Che  average  of  about  0.1  Ry. 
On  the  other  hand,  all  the  a^^  states  have  been  affected  irregularly  by 
the  Pb  interstitial.  In  order  to  describe  their  characteristics,  we 
show  their  partial-wave  coefficients  and  electronic  charge  distributions 
in  Table  III  (for  comparison,  the  5C2  bottom  conduction  band  state  is 
also  included) . 

The  ISj^  state  represents  a bonding  combination  of  the  Te  5s,  IPb 
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(Interstitial  Fb)  6s,  and  the  Pb  6s  atomic  basis  states  vhlle  the  2a^ 
state  represents  an  antibonding  combination  of  the  Te  5s-p  hybrid  and 
the  IFb  6s  basis  states.  The  3a^  state  corresponds  to  the  2a^  state  of 
the  OPb^Te^  cluster;  however,  its  eigen-energy  has  not  been  lowered  as 
much  as  the  e,  t^,  or  states  because  the  IFb  6s  basis  state  is  combined 
with  the  Pb  6s-p  hybrid  state  to  produce  an  antibonding  behavior  for 
the  3a^  state. 

The  4aj^  state  exhibits.  Interestingly,  a rather  strong  antibonding 
between  the  Te  5p  and  the  IPb  6s  basis  states;  thus,  its  level  has 
been  raised  from  the  corresponding  3a^  level  in  the  top  valence  band 
of  the  OPb^Te^  cluster  to  the  conduction  band  of  the  present  Pb 
interstitial  cluster.  The  contSur  map  of  its  wave  function  in  Fig.  5 
clearly  shows  such  antibonding  character.  The  level  is  higher  than 
the  5t2  level  the  electron  state  of  which  corresponds  to  the  5t2  state 
of  the  OPb^Te^  cluster  and  hence  represents  the  bottom  of  the  conduction 
band. 

As  we  introduce  the  Pb  interstitial  in  the  OPb^Te^  cluster,  we  get 
an  extra  a^  state  in  the  combination  of  the  la^  and  2a.^  states  and  lose 
an  state,  which  corresponds  to  the  4aj^  state  of  the  OPb^Te^  cluster, 
from  the  top  valence  band  to  the  conduction  band.  The  net  result  is 
that  we  have  the  same  number  of  electron  levels  in  the  valence  bands  of 
the  Pb  interstitial  cluster  as  in  the  case  of  the  OPb/Te,  cluster. 

However,  in  the  Pb  interstitial  cluster,  we  have  four  more  valence  electrons 
due  to  the  interstitial  Pb  atom.  Therefore,  we  are  left  with  four  extra 
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electrons  In  the  bottom  conduction  band  state  St^*  Thus  a Pb  inter- 
stitial gives  four  free  electrons  to  the  PbTe  crystal. 

The  constant  potential  of  the  Interatomic  region,  near  which  we 
suppose  the  crystal  bandgap  to  be  located.  Is  only  slightly  below  the 
St2  level.  However,  the  gap  between  the  4t2  and  the  5t2  levels  Is 
large  enough  to  assure  that  the  bandgap  of  the  crystal  Is  Indeed  located 
in  the  gap  and  consequently  the  5t2  level  beongs  to  the  conduction  band. 

Commonly  It  has  been  thought  that  a Pb  Interstitial  would  give  up 

2+ 

only  two  free  electrons  by  becoming  a Pb  ion  and  depositing  two 
electrons  in  the  conduction  band  of  the  crystal.  Such  a model  was 
made  apparently  without  considering  the  antibonding  effect  due  to  the 
interstitial  which  drives  up  a valence  band  state  of  the  perfect  crystal 
to  the  conduction  band. 

Now  let  us  look  into  the  Te  interstitial  problem.  The  la^^  and  la.^ 

states  of  the  Te  interstitial  cluster  originate  from  the  la^^  state  of 

the  OPb.Te,  cluster  and  ITe  (interstitial  Te)  5s  atomic  orbital.  The 
4 4 

3a^  state  is  the  perturbed  state  originating  from  the  Za^^  state  of 
the  OPb^Te^  cluster.  The  4a^  state  corresponds  to  the  Sa^^  state  of  the 
OPb^Te^  cluster,  but  due  to  its  antibonding  character  with  the  ITe  5s 
basis  state  its  eigen-energy  has  been  raised  unlike  the  other  a^  states 
Each  of  the  e and  t^^  states  of  the  Te  interstitial  cluster  cor- 
responds to  the  state  with  the  same  denomination  in  the  OPb^Te^  cluster. 
So  does  the  lt2  and  2t2  states.  However,  the  situation  is  not  very 
clear  for  the  3-6  t2  states. 


16 


Th«  state  seems  to  originate  from  the  3t2  state  of  the  OPb^Te^ 

cluster,  but  It  has  a significant  component  of  the  Pb  6s  basis  state  as 

Che  2t.  state  of  the  OPb.Te,  cluster  has.  The  4c.  and  the  6c.  states 
2 4 4 2 2 

seem  to  correspond  to  the  4C2  and  bc^  states  of  the  OPb^Te^  cluster 
respectively.  But  the  Sc^  state  does  not  correspond  to  any  one  particular 
state  of  the  OPb^Te^  cluster;  it  seems  to  be  a perturbed  com- 

bination of  all  Che  c^  states  in  Che  cop  valence  band  and  the  conduction 
band.  It  is,  therefore,  not  certain  where  the  crystal  bandgap  would  be 
located  in  relation  to  the  Sc^  level. 

If  we  assume  chat  the  bandgap  of  the  crystal  is  still  very  near  the 
constant  potential  in  the  interatomic  region,  the  St^  level  belongs  to 
Che  valence  baud.  Then  a Te  interstitial  would  give  two  free  holes  to 
Che  crystal  since  the  Sc^  level  is  occupied  by  only  four  electrons.  This 
picture  implies  chat  all  the  5s  and  5p  atomic  levels  of  the  Te  inter- 
stitial bd.arg  to  Che  cop  valence  band.  If  the  Te  interstitial  ever 
exists  in  PbTe,  the  interstitial  complex  will  be  greatly  relaxed  due  to 
the  large  size  of  the  Te  ion;  this  fact  further  obscures  the  conclusion 
drawn  here. 


N 
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V.  Pb  AND  Te  FRENKEL  CEFECTS 

A Freakel  defect  in  a crystal  Is  formed  when  a coastltuent  atom  Is 
displaced  from  its  regular  site  In  the  lattice  to  an  Interstitial  site. 
Therefore,  It  may  be  considered  as  an  aggragate  of  a vacancy  and  an 
Interstitial  of  the  same  atomic  species.  Such  an  aggragate  leads  to  many 

different  types  of  Frenkel  defect  complexes  depending  on  the  separation  of 
the  interstitial  and  the  vacancy.  In  this  section  we  will  deal  only 
with  a Frenkel  defect  in  which  the  pair  are  placed  as  close  as  possible. 

The  Frenkel  defect,  then,  can  be  conveniently  modeled  within  the  OPb^Te^ 
cluster  by  displacing  an  atom  to  the  center  of  the  cluster.  This  de- 
creases the  symmetry  of  the  complex  from  the  T^  to  the  point  group 
and  reproduces  the  correct  symmetry  of  such  a defect  in  a crystal. 

The  Pb  Frenkel  defect  is  another  kind  of  important  point  defect 
species  in  PbTe.  But,  the  Te  Frenkel  defect  seems  to  be  difficult  to 
be  formed  for  the  same  reason  given  to  the  Te  interstitial. 

The  electron  energy  levels  of  the  Pb  and  Te  Frenkel  defect  clusters 
are  shown  in  Fig.  6,  compared  with  those  of  the  OPb^Te^  cluster. 

We  will  discuss  the  Pb  Frenkel  defect  first.  The  la^^  and 
states  are  the  bonding  and  antibonding  combinations,  repectively,  of  the 
Te  Ss  and  the  IPb  6s  basis  states.  Between  the  two  levels  there  are  two 
closely  lying  levels,  2a and  le,  which  originally  belonged  to  the  Te  s 

band  of  the  OPb.Te,  cluster.  The  2e  and  Aa.  levels  now  constitute 
A A 1 

the  perturbed  Pb  s band.  Next,  the  levels  from  5a^  to 


6e  make  up 
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the  cop  valence  band,  and  the  levels  of  7e  and  up  the  conduction  band. 

One  thing  to  be  noted  here  Is  that  the  8a^  state  is  an  antibonding 
combination  of  the  Te  Sp  and  IFb  6s  basis  states,  thus  exhibiting  similar 
characteristics  as  the  4a^  state  of  the  Pb  Interstitial  cluster.  The 
8a^  level  In  the  conduction  band  Is  the  only  level  vhlch  has  been  re- 
moved from  the  valence  bands  uhlle  no  electrons  are  lost  as  the  Pb 
Frenkel  defect  is  formed.  Therefore,  the  bottom  conduction  band  state 
7e  carries  tvo  electrons  which  will  become  free  in  the  crystal. 

The  result  of  the  Pb  Frenkel  defect  is  exactly  an  algebraic  sun 
of  the  results  of  the  separate  Pb  vacancy  and  Pb  laterstital.  In  the 
Pb  s band,  the  missing  Pb  6s  level  due  to  the  Pb  vacancy  is  compensated 
by  the  6s  level  of  the  Pb  Interstitial,  but  the  Pb  interstitial  in  turn 
loses  an  a^  state  from  the  top  valence  band  to  the  conduction  band. 
Therefore,  we  may  say  that  the  two  free  holes  of  the  Pb  vacancy  are 
cancelled  by  the  two  out  of  the  four  free  electrons  of  the  Pb  inter- 
stitial, thereby  the  Pb  Frenkel  defect  leaves  two  free  electrons  for  the 
crystal. 

Now  let  us  turn  to  the  Te  Frenkel  defect.  The  la^^,  le,  and  2a 
levels  constitute  the  perturbed  Te  s band,  and  the  levels  from  3a^  to 
4a^  constitute  the  perturbed  Pb  s band.  Next,  the  set  of  levels  from 
3e  to  makes  up  the  perturbed  top  valence  band  while  the  levels  Sa^^ 
and  up  correspond  to  the  conduction  band  of  the  Te  Frenkel  defect  cluster. 
Therefore,  the  bandgap  in  the  cluster  model  is  between  the  7a^  and  8a^ 
levels.  In  this  Te  Frenkel  defect  cluster,  the  constant  potential  in  the 
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lateracomlc  region  does  not  correctly  Indicate  the  location  of  the 
crystal  bandgap. 

There  is  ambiguity  in  Interpreting  the  6e  level  which  Is  a perturbed 
combination  of  some  conduction  band  states  and  some  valence  band  states. 
However,  let  us  assume  that  it  belongs  to  the  valence  band  since  Its 
eigen-energy  Is  equal  to  that  of  the  top  valence  level  This  makes 

all  the  valence  band  levels  filled  and  all  the  conduction  band  levels 
empty;  consequently  the  Te  Frenkel  defect  does  not  contribute  any  to  the 
conductivity  of  the  PbTe  crystal.  This  result  Is  also  an  algebraic 
sum  of  the  separate  results  from  a Te  vacancy  and  a Te  Interstitial. 

The  Frenkel  defect  of  Pb  or  Te  which  we  have  illustrated  thus  far 
has  the  vacancy  and  the  interstitial  as  close  to  each  other  as  possible. 
Tet  they  shew  that  the  result  Is  the  same  as  the  sum  of  the  individual 
effects  from  the  vacancy  and  the  interstitial.  Therefore,  any  Frenkel 
defect  regardless  of  the  separation  between  the  pair  will  give  the  same 
result  because  the  interaction  between  the  vacancy  and  the  Interstitial 
becomes  weaker  as  their  separation  gets  larger. 
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VI.  CONCLUSION  AND  FURTHER  DISCUSSION 

Wc  suxDmarlze,  In  the  cable  below,  the  type  and  Che  number  of  free 
carriers  chat  each  point  defect  species  would  Induce  in  the  PbTe  crystal. 

Effects  of  the  point  defects  on  the  conductivity  of  PbTe,  as 
deduced  from  the  SCF  Xa  cluster  model  calculations. 


Defects  Species  Pb 


Te 


Vacancy 
Interstitial 
Frenkel  defect 


2 holes  2 electrons 

4 electrons  2 holes  (?) 

2 electrons  No  carrier  (?) 


The  difficulties  we  have  encountered  in  interpreting  the  energy 
levels  for  the  Te  interstitial  and  the  Te  Frenkel  defect  may  have  some 
relationship  with  the  unlikeliness  of  their  occurrence  in  the  PbTe  crystal. 
If  we  disregard  the  two  point  defect  species,  a p-type  PbTe  semiconductor 
can  be  made  only  through  the  Pb  vacancies,  as  all  the  other  kinds  lead 
to  the  n-type  semiconductor. 

The  results  we  obtained  for  the  Pb  and  Te  vacancies  in  the  PbTe 
crystal  were  shown  to  agree  with  the  results  from  perturbative  treatments. 
However,  the  conclusions  for  the  interstitials  and  the  Frenkel  defects 
arc  not  in  agreement  with  such  strictly  perturbative  considerations  made 
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by  one  of  the  authors. In  the  present  cluster  model  calculations, 
ve  showed  chat  there  are  some  electron  levels  which  are  affected  by 
Che  perturbing  potential  in  the  opposite  direction  due  to  particular 
anclbonding  mechanisms  although  most  energy  levels  follow  the  perturbing 
potential. 

The  defect  clusters  for  the  Pb  vacancy,  both  kinds  of  interstitials, 
and  the  Pb  Frenkel  defect  are  open-shell  systems.  However,  spin- 
polarized  calculations  have  not  been  pursued,  because  the  incompletely 
filled  levels  do  not  represent  localized  states  in  the  bandgap  of  Che 
crystal,  but  they  fall  into  one  energy  band  or  another.  The  electrons 
or  holes  in  the  energy  bands  will  become  free  in  the  crystal. 

The  effects  of  the  intrinsic  point  defects  in  PbTe  are  very  much 
localized  in  space.  Therefore,  such  small  clusters  as  have  been  used 
in  the  present  paper  were  good  enough  to  show  the  relevent  short-range 
effects  of  the  defects.  In  fact,  they  lead  to  Che  same  results  as 
obtained  from  larger  cluster  models  by  Hamstreet^^  in  the  case  of 
vacancies . 
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APPENDIX:  MASS  VELOCITY  AND  DARWIN  RELATIVISTIC  CORRECTIONS  IN  THE  ATOMIC 
RADIAL  WAVE  EQUATION 


The  two-coapoaenc  relativistic  Pauli  equation  for  an  electron  Is 
obtained  by  eliminating  tvo  small  components  from  the  four-component  Dirac 
equation.  Specializing  to  an  electron  in  a spherically  symmetric  potential, 
V(r),  one  gets  the  radial  Pauli  equation, 

((i  ♦ T (-  7 ^ ^ * '''■■> 

- “x  i ^ ^ * [r:  i])l 


- E (R;E), 


where  a is  the  fine-structure  constant,  (137.037)  in  a.u.,  and  E is  the 

_2 

eigen-energy  of  an  electron  less  its  rest  energy,  2a  a.u.  In  the  above 

Eq.  (1)  , the  upper  quantity  (-1)  refers  to  the  case  when  the  total  angular 

momentum  J - 2.  + 1/2,  and  the  lower  quantity  (2,  + 1)  to  the  case  when 

J ■ 2,  - 1/2  for  an  electron  in  the  atomic  orbital  of  angular  momentim  2. 

Rather  than  solving  Eq.  (1)  directly,  we  will  proceed  to  get  a simplified 

2 

equation  which  retains  only  the  terms  of  the  order  or  unity  and  a . This 
equation  can  be  written  as 

tH^(r)  + H (r)  + H, (r)  + H_^(r)]  R,  (r;E)  - E R,  (r;E).  (2) 

o m a so  jc  x. 

In  the  above  equation,  H (r)  represents  the  non-relatlvistic  Hamiltonian, 

o 
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\(r)  the  relativistic  mass-velocity  term,  (r)  the  relativistic  Darwin 

term,  and  H (r)  the  spin-orbit  coupling  term;  they  are  expressed  by 
so 


H (r) 
o 

H (r) 
m 

H^(r) 


.•-L  A ^2  + v(r) 

2 dr  ^ dr  ^ ^ 2 

r r 

«2  2 
[E  - V(r)]^  , 

_ 2I  d^(r)  _d_ 

4 dr  dr  ’ 


(3) 

(4) 

(5) 


and 


so 


4 

0 


' -2  *1  i dV(r) 

.1  + ij ' 


for  2^1, 
for  2-0. 


(6) 


V. 

In  crystals  and  molecular  clusters,  electron  states  are  classified 

by  representations  of  a symmetry  group.  The  inclusion  of  spin-orbit  coupling 

requires  the  use  of  double  group  representations  while  the  SCF  Xa  multiple- 

scattering  method  is  based  on  the  symmetry  of  single  groups.  Therefore, 

in  the  present  modification  for  relativistic  corrections,  we  drop  the  spin- 

orbit  coupling  term.  The  fully  relativistic  formalism  for  the  multiple- 

25 

scattering  method  has  been  recently  worked  out  by  Yang  and  Rabii. 

Since  we  use  the  muffin-tin  approximation  for  tlje  multiple-scattering 
method,  the  relativistic  mass-velocity  and  Darwin  terms  can  be  incorporated 
as  they  are  shovm  in  Eqs.  (4)  and  (5).  These  terms  are  necessary  only 
for  heavy  atoms.  In  implementing  such  corrections, 

we  will  directly  solve  Eq.  (2)  less  the  spin-orbit  coupling  term.  Before 
we  describe  the  procedure  in  solving  the  equation,  let  us  briefly  review 


the  non-relativistic  case. 
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The  non-relaciviscic  radial  wave  equation  is 


R^(r;E)  - E Rj^  (r;E). 


Nov  by  the  use  of 


(r;E)  - r Rj^  (r;E), 


Eq.  (7)  can  be  transforxed  into 

.2 


JI.VM 

dr  r 


Pj^(r;S)  - 0 


which  can  be  rewritten  as 


(7) 


(8) 


(9) 


where 


d (r;E) 
dr 


gj^Cr-.E) 


(r;E)  (r;E), 


E - V(r)  - 


l(i  + 1) 


(10) 


(11) 


Equation  (10)  can  be  solved  numerically  by  the  Numerov  method  with  a r- 

mesh,  given  the  two  initial  values  (r^;E)  and  ?j^(r2;E).  The  Nunerov 

method  and  the  scheme  of  starting  numberical  integration  is  described 

26 

by  Mattheiss,  et  al. 

Now  adding  the  relativistic  mass-velocity  and  Darwin  terms,  the  radial 
wave  equation  becomes 

f;7  37r2^+V(r)  [E  -V(r)l  = 


(12) 


-T^  ^\h 


Let  us  define  chat 


(r;E)  - V(r)  -t-  ^ J [r  - V(r)  ] “ - - i-E  . 


V-.-/  -vw  ^ _2  4 r“  r dr 


then  Eq.  (13)  becomes 


U' 


;E)j  (r;! 


:)  - 0. 


where 


f(r)  - ^ V'(r) 


Substituting 


(r;E)  - y(r)  0^  (r;E), 


In  Eq.  (15),  we  now  have 


-1  - 2(r-E) 
7(r)  dr 


^ y”(r)  -^  f(r)  v'(r)  ^ q 

y(r) 


(14) 


(15) 


(16) 


(17) 


(18) 


In  order  to  eliainate  the  second  tern  in  the  above  equation,  we  choose  the 


arbitrary  function  y(r)  such  that 


4 
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2y'(r)  + f(r)  y(r)  - 0. 


aasiming  chat  y(r)  does  not  have  any  zeros.  The  solution  of  Eq.  (19)  is 


(19) 


y(r)  - exp 


piv<o]. 


(20) 


which  indeed  does  not  vanish  for  any  r. 

Equation  (18)  can  now  be  written  as 


where 


Qj^"(r;E)  - (r;E)  (r;E), 


G (r-E)  - g (r-E)  - 

2 ^ Id  + 1) 


(21) 


- (1  + ^ E)V(r)  - ^ [V(r)]"  + 

r*" 

2 2 . 4 

-E  + ^ V"(r)  + ^ p V'(r)  + |^  [V'(r)]^. 


After  evaluating  G (r;E)  of  Eq.  (22)  and  solving  Eq.  (21)  for 
Q (r;E)  by  the  Nunerov  tnethod  as  in  the  non-relativistic  case,  Che  radial 

X( 

wave  function  R (r;E)  can  be  found  by 
X* 


(22) 


(r;E)  - ^ exp^- V(r)j  (r;E). 


(23) 


It  was  found  that  an  iterative  scheae  is  very  efficient  in 
obtaining  the  two  initial  values,  (rj^;E)  and  (r2;E),  which  are 
necessary  to  start  the  numberical  integration  for  (r;E). 

The  scheae  is  explained  as  follows.  For  the  several  starting 


mesh  points  for  r 0,  V(r)  -TZlr,  Therefore,  in  the  expression  for 

(r;E)  in  Eq.  (22),  the  terms  involving  E are  insignificant  at  r'XD  even  for 
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‘ core  electron  levels.  We  now  set  (r2)  ® (r'2;E)  ■ 1.  since  the 

radial  wave  function  (r;E)  can  be  normalized  later.  First,  assuming 

(rj^)  ■ 1,  we  use  the  Numerov  method  to  compute  (r)  for  r ■ r^,  r^,  and 
rj.  For  these  mesh  points,  R^  (r)  ■ Rj^  (r;E)  can  be  obtained  Immediately 
by  Eq.  (23).  Extrapolation  from  the  fixed  R^  (r2)  and  R^^  (r)  for  r ■ r^. 
and  r^,  we  get  Rj^  (rj^)  and  subsequently  (r^^)  which  is  now  compared  with 

the  value  with  which  we  started  the  cycle.  The  cycle  is  Iterated 
until  the  two  consecutive  values  for  (r^^)  are  close  enough  to  make  the 
relative  error  smaller  than  a pre-set  tolerance.  The  relative  error  can 
be  smaller  than  10  ^ very  quickly  only  after  a few  iterations. 

In  this  Iterative  scheme  for  the  starting  values,  we  used 

' Rj^  (r)  Instead  of  (r)  for  the  extrapolation  step,  because  the  former  is 
a more  smooth  function  than  the  latter. 

As  the  starting  values  are  Independent  of  E,  we  need  to 

calculate  them  only  once  for  a given  potential  energy  in  searching  for  the 
electron  eigen-energies. 
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TABLE  I 


Electron  states  of  the  OPb^Te^  cluster. 

They  are  listed  in  the  order  of  increasing  eigen-energy.  Only  the  important 
components  of  the  normalized  partial-wave  coefficients  and  the  electronic 
charge  distributions  are  shown.  (The  "In"  and  "Out"  represent  the  inter- 
atomic and  the  extramolecular  region  of  the  cluster  respectively.)  The  con- 
stant potential  in  the  interatomic  region  is  -0.335  Ry. 


States 

Eigen- 
energy  (Ry) 

Partial-wave 

coefficients 

Electronic  Charge 
distribution 

^*1 

-1.095 

Te  s 
0.89 

Pb  s 
-0.35 

Te 

0.80 

Pb 

0.15 

-1.050 

Te  s 
0.95 

Pb  s 
0.14 

Te 

0.91 

Pb 

0.04 

-0.766 

Pb  s 
0.87 

Te  p 
-0.30 

Pb 

0.77 

Te 

0.15 

-0.737 

Pb  s 
0.78 

Te  s 
0.40 

Te  p 
0.39 

Pb 

0.61 

Te 

0.31 

le 

-0.519 

Te  p 
0.77 

Pb  p 
0.48 

Te 

0.60 

Pb 

0.23 

In 

0.16 

-0.502 

Te  p 
-0.63 

Pb  p 
-0.50 

Pb  s 
0.38 

Te 

0.40 

Pb 

0.40 

In 

0.09 

-0.462 

Te  p 
0.83 

Pb  p 
-0.42 

Te 

0.69 

Pb 

0.17 

In 

0.12 

-0.450 

Te  p 
0.72 

Pb  p 
-0.40 

Te  p 
0.38 

Te 

0.65 

Pb 

0.20 

In 

0.14 

^^2 

-0.401 

Te  p 
0.77 

Pb  p 
0.29 

Te 

0.68 

Pb 

0.20 

In 

0.08 

-0.211 

Pb  p 
0.70 

Out  p 
0.27  ■ 

Te  p 
-0.20 

Pb 

0.51 

In 

0.32 

Te 

0.08 

Out 

0.08 

'h* 

-0.155 

Pb  p 
0.73 

Te  p 
0.39 

Pb 

0.53 

In 

0.30 

Te 

0.15 

6t^* 

-0.147 

Pb  p 
0.65 

Telp 

0.38 

Out  d 
0.36 

Pb 

0.50 

Te 

0.18 

In 

0.17 

Out 

0.14 

2c  * 

-0.106 

Pb  p 
0.66 

Te  p 
-0.51 

Pb 

0.43 

Te 

0.26 

In 

0.28 

* 

-0.098 

Pb  p 
0.60 

Te  p 
-0.51 

Out  f 
0.38 

Pb 

0.39 

Te 

0.34 

Out 

0.15 

In 

0.09 

Excited  empty  electron  states. 
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TABLE  II 


D«scripcion  of  the  6aj^  and  6e  states  of  the  Te  vacancy  cluster. 

Only  the  Inportant  components  of  the  normalized  partial-wave  components 
and  the  electronic  charge  distributions  are  shown.  (The  "In", 

and  "Out"  represent  the  Te  vacancy  center,  ineteratomic , and  extra- 
molecular  regions  of  the  cluster  respectively.) 


States 

Eigen- 
energy  (Ry) 

Partial-wave 

coefficients 

Electronic  charge 
distribution 

6ai 

-0.323 

Pb  p 

Te 

Te  p 

P 

Te 

Pb 

In 

^Te 

0.70 

-0.39 

-0.22 

0.13 

0.50 

0.23 

0.17 

6e 

-0.245 

Pb  p 

P 

Te 

Te  p 

Cut  d 

Pb 

In 

^Te 

0.77 

-0.23 

0.16 

-0.09 

0.60 

0.28 

0.05 

« 
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TABLE  III 

Descrlpcloo  of  che  l-4a^  and  Sc^  states  of 
the  Pb  interstitial  cluster. 


Only  the  Important  components  of  the  normalized  partial-wave 
coefficients  and  electronic  charge  distributions  are  shown. 
(The  "IPb”  and  "In"  represent  the  Pb  interstitial  center  and 
che  interatomic  region  of  the  cluster  respectively.)  The  con- 
stant potential  in  the  interatomic  region  of  the  cluster  is 
-0.351  Ry. 


States  Eigen- 

. engergy  (Ry) 

Partial-wave 

Coefficients 

Electronic  charge 
d istribution 

^*1 

-1.361 

Te  s 

IPb  s 

Pb  s 

Te 

IPb 

Pb 

1 

0.66 

-0.53 

-0.40 

0.50 

0.26 

0.24 

-1.013 

Te  s 

IPb  S 

Te  p 

Te 

IPb 

X 

0.72 

0.45 

-0.39 

0.67 

0.20 

3*1 

-0.777 

Pb  s 

Pb  p 

IPb  s 

Pb 

In 

IPb 

X 

0.83 

0.31 

-0.29 

0.79 

0.09 

0.09 

4*1 

-0.333 

Te  p 

IPb  s 

Te 

IPb 

In 

X 

0.79 

0.40 

0.66 

0.16 

0.14 

-0.345 

Pb  p 

Te  p 

Pb 

In 

Te 

4 

0.74 

0.28 

0.56 

0.23 

0.13 

Fig.  1 Cluster  model  of  OPb^Te^  for  the  PbTe  crystal. 

Fig.  2 One-electroQ  energy  levels  of  the  OPb^Te^  cluster  (a),  and 
the  energy  bands  of  the  PbTe  crystal  (b)  as  calculated  by 

Lin  and  Klelnman  (ref.  21).  — represents  the  constant 

potential  In  the  interatomic  region  of  the  cluster.  The 
crystal  energy  b^nds  were  shifted  so  that  the  centers  of  the 
band  gaps  match  In  the  two  energy  level  structures. 

Fig.  3 One-electron  energy  levels  for  the  Pb  vacancy  cluster  (a) 

and  the  Te  vacancy  cluster  (c) , compared  with  the  levels  for 

the  OPb^Te^  cluster  (b) . represents  the  constant 

potential  in  the  Interatomic  region  of  a cluster,  and 
indicates  the  highest  occupied  level  with  its  occupation 
number  in  the  parenthesis. 

Fig.  4 One-electron  energy  levels  for  the  Pb  interstitial  cluster 
(a)  and  the  Te  interstitial  cluster  (b) , compared  with  the 
levels  for  the  OPb^Te^  cluster  (c) . — • — • — represents 

the  constant  potential  in  the  interatomic  region  of  a cluster, 
and  indicates  the  highest  occupied  level  with  its 

occupation  numbuer  in  the  parenthesis. 

Fig.  S Contour  map  of  the  4a^  xrave  function  of  the  Pb  interstitial 
cluster.  It  was  plotted  in  a plane  which  bisects  the  two  Te 
atoms  at  the  right,  the  two  Pb  atoms  at  the  left,  and  the 


Interstitial  Pb  atom  at  the  center.  The  broken-line  circles 


are  the  projections  of  Che  spheres  of  the  atomic  centers 
and  Che  outersphere  of  the  cluster.  The  contour  curves 
were  selected  such  Chat  the  logarithms  of  the  magnitudes  of 
the  wave  function  values  are  equally  spaced. 

Fig.  6 One-electron  energy  levels  for  the  Pb  Frenkel  defect  cluster 
(a)  and  the  Te  Frenkel  defect  cluster  (b) , compared  with 

the  levels  for  the  OPb^Te^  cluster  (c).  respresents 

the  constant  potential  in  the  interatomic  region  of  a cluster, 
and  'V'WN  indicates  the  highest  occupied  level  with  its 
occupation  number  in  the  parenthesis. 
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ABSTRACT 


A GaP  crystal  was  simulated  by  the  smallest-possible  neutral  cluster 
model,  GaPHg,  which  is  comprised  of  a GaP  molecule  and  six  boundary  hydrogen 
atoms  replacing  three  gallium  and  three  phosphorus  nearest-neighbor  atoms. 
The  results  of  the  calculation  by  the  SCF  Xs  multiple-scattering  method 
of  Johnson  showed  that  the  energy-level  structure  of  the  cluster  fairly 
well  duplicates  the  energy-band  structure  of  the  GaP  crystal.  In  this 
cluster  model  some  hydrogenic  surface  electron  states  were  found  and  they 
were  made  unoccupied  in  the  SCF  iterations.  Next,  the  phosphorus  atom  of 
the  cluster  was  substituted  by  a nitrogen  atom,  resulting  in  a cluster  of 
GaNHg,  to  simulate  the  GaP  crystal  embedded  with  a nitrogen  impurity.  Care- 
ful analysis  of  the  electronic  structure  of  the  GaNH,  cluster  in  comparison 

0 

with  that  of  the  GaPH^  cluster  reveals  a bound  electron  state  in  the  bandgap 

o 

which  represents  the  isoelectronic  electron  trap  due  to  the  nitrogen  impurity 
In  the  Gap  crystal.  The  present  cluster  model  calculations  showed  that 
the  bound  electron  state  originates  from  the  state  of  the  conduction 
bend,  thus  explaining  the  occurrence  of  the  efficient  luminescence  A line 
due  to  isoelectronic  nitrogen  impurities  in  GaP. 
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I.  INTRODUCTION 

The  Group  III-V  semiconductors  have  been  studied  extensively  over 
the  last  decade  due  to  their  Important  technological  Implications  In  electronic 
devices.^  In  most  of  the  applications,  the  materials  are  doped  by  impurities 
or  by  Intrinsic  defects  in  order  to  control  the  electrical  and  optical 
properties. 

The  shallow  hydrogenlc  impurity  electron  states  in  semiconductors  in  general 

are  now  well  understood,  both  experimentally  and  theoretically.  However,  for 

the  localized  bound  electron  states  caused  by  some  point  defects  in  semiconductors 

there  are  still  a lot  more  to  be  comprehended.  In  this  area  theoretical 

knowledge  is  particularly  lacking;  there  have  been  only  several  first- 

principles  calculations  for  the  kind  of  point  defects,  such  as  the  vacancies 
2 3 

in  silicon  and  lead  tellur ide  by  the  energy-band  method  through  the  Green's 

4 

function  formalism  of  Roster  and  Slater  , and  the  nitrogen  impurity  and 
vacancy  in  diamond^  by  the  LCAO-MO  method.  Recently  the  SC?  Xa  multiple- 

scattering  method  of  Johnson^  was  applied  to  the  vacancies  in  lead  chalco- 

7 8 

genldes  by  Hamstreet  and  to  the  impurities  in  silicon  by  Cartllng  to  find 

the  electronic  structures  of  the  point  defect  complexes  by  the  use  of  cluster 

models.  This  method  is  very  efficient,  but  careful  consideration  must  be 

given  to  the  configuration  and  the  boundary  condition  of  a cluster  model  in 

order  to  simulate  a crystal  successfully.  The  authors  applied  the  method, 

a 

In  a previous  paper,  to  simulate  a LIF  crystal  by  a smallest-possible  neutral 
cluster  model  and  obtained  the  bulk  properties  of  the  crystal  in  very  good 
agreement  with  experiments.  The  F-center  problem  of  the  LIF  crystal  was  also 
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studied.  Thea  In  a subsequent  paper,  the  same  formalism  was  used  for  a 
PbTe  crystal  to  find  the  effects  of  various  intrinsic  point  defects  to  the 
conductivity  of  the  PbTe  crystal. 

In  this  paper,  we  extended  our  study  to  GaP  crystal,  which  is  one  of 

the  Group  IIl-V  semiconductors.  The  crystal  was  first  simulated  by  the 

smallest-possible  neutral  cluster  model,  and  then  the  isoelectronic  electron 

trap  caused  by  a nitrogen  impurity  was  investigated  with  the  same  cluster 

format.  The  nitrogen  isoelectronic  electron  trap,  which  was  first  discovered 

by  Thomas,  et  al.^^  is  responsible  for  the  efficient  green-light  luminescence 

in  the  CaP  light-emitting  diode.  This  problem,  as  far  as  we  are  aware  of, 

has  not  been  treated  by  any  theoretical  means.  Since  the  materials  we 
• 9 10 

studied  in  the  previous  papers  ’ are  very  much  ionic,  the  highly  covalent 
GaP  here  will  provide  a complementary  example  in  the  treatment  of  crystals 
and  their  point  defects  by  the  SCF  Xx  multiple-scattering  method. 

The  GaPH,  cluster  model  for  a GaP  crystal  is  created  in  Section  II,  and 

Q 

the  effects  of  the  nitrogen  impurity  in  the  cluster  is  explained  in  Section 
III.  Finally,  in  Section  17  a summary  is  given. 


II.  GaP  CRYSTAL 

The  GaP  crystal  has  the  zinc-blende  structure,  whose  symmetry  belongs 

2 *12 
CO  the  space  group  td  , with  a nearest-neighbor  distance  of  2.36  A.  It 

has  an  indirect  bandgap  of  0.172  Ry  (2.34  eV)  at  0*K  which  decreases  to 

0.166  Ry  (2.26  eV)  at  300*K.^ 

In  a simple  valence-bond  picture  of  GaP,  each  atom  is  described 
as  having  four  covalent  bonds  tetrahedrally  directed  to  its  four  nearest 
neighbors.  Once  electron  of  each  P atom  is  transferred  to  each  Ga  atom  in 
this  model  so  that  every  atom  may  have  four  equivalent  hyrid  orbitals  made 
of  atomic  a and  o functions. 

A.  Cluster  Model  for  the  GaP  Crystal  • ^ ' 

If  we  choose  the  smallest-possible  cluster  model  for  the  GaP  crystal,  | 

i 

retaining  the  full  point  group  symmetry  of  the  crystal  while  having  a 
phosphorus  atom  at  the  center  of  the  cluster  so  that  it  provides  a con- 
venient site  for  a nitrogen  Impurity,  it  will  be  a ?Ga^  cluster  where  four 
gallium  atoms  are  placed  tetrahedrally  around  the  phosphorus  atom.  This 
cluster  is  very  much  unbalanced  regarding  the  composition  of  the  Ga?  crystal; 

since  the  cluster  contains  only  one  phosphorus  atom  with  four  gallium  atoms, 

14 

the  energy-bands  of  the  cluster  will  have  less  contribution  from  the 
phosphorus  chan  the  gallium  atoms  compared  to  the  energy-bands  of  the  crystal. 

Also,  Che  electronic  charge  transfer  between  the  phosphorus  and  gallium 
atoms  in  the  cluster  will  not  render  the  same  meaning  as  that  in  the  crystal. 

A SCF  ^ mulclple-scatcerlng  calculation  for  the  PGa^  cluster  model 
shows  chat  there  is  one  electron  left  over  after  Che  valence  bands  are 
completely  filled.  This  is  explained  as  follows.  The  top  valence  band  level 
of  Che  cluster  is  a six-fold  degenerate  t2  level  whose  wave  function  is 
composed  primarily  of  P 3p. atomic  orbitals.  Since  there  are  only  three 
P 3p  electrons,  the  rest  of  the  C^  level  oust,  be  filled  by  three  Ga  4p 
electrons.  However,  there  are  four  Ga  4p  electrons  in  the  cluster,  one  for 


each  Ca  atom,  and  the  one  remaining  electron  goes  over  to  a conduction  band 

level  whose  wave  function  Is  made  primarily  of  Ga  4p  atomic  orbitals. 

Therefore,  we  have  to  remove  this  conduction  band  electron  so  that  ve  have 

completely  filled  valence  bands  and  completely  empty  conduction  bands,  as 

In  the  case  for  the  perfect  GaP  crystal.  This  means  that  ve  have  to  freeze 

the  conduction  band  electron  on  a Watson  sphere, thus  resulting  in  a 

(PGa^)'*’  cluster.  The  use  of  the  Watson  sphere  will  make  the  cluster  model 

more  artifical.  Furthermore,  the  calculation  of  the  (PGa,  )"*"  cluster  shows 

4 

no  aj^  state  in  the  conduction  band  which  would  correspond  to  the  state 
of  the  GaP  crystal,  the  buttom  conduction  band  state  at  the  center  of  the 
Brillouin  zone.  Therefore,  no  detailed  comparison  of  the  band-edge  structure 
can  be  made  between  the  cluster  and  the  crystal,  and  this  makes  further 
investigation  of  the  isoelectronic  electron  trap  meaningless. 

In  order  to  accommodate  all  the  electrons  equally  without  the  need 
of  freezing  any  electrons,  the  cluster  model  must  be  neutral,  i.e.  it  must 
contain  the  constituent  atoms  of  the  crystal  in  the  same  proportion  as  the 
crystal.  The  smallest  neutral  cluster  for  the  GaP  crystal,  but  larger  than 
a mere  molecule,  is  a Ga^P^  cluster  in  which  each  atom  of  an  inner 
GaP  molecule  has  four  nearest-neighbor  atoms  of  the  other  kind.  The  cluster 
symmetry  now  belongs  to  the  point  group,  a subgroup  of  the  Td  group  of 
the  crystal.  It  is  impossible  to  formulate  a cluster  which  is  both  neutral 
and  fully  symmetric.  For  the  cluster,  however,  there  will  be  a lot 

of  dangling  bonds  from  the  six  outer  atoms  due  to  the  covalent-bonding 
nature.  In  fact,  a SCF  Xa  multiple-scattering  calculation  for  the  Ga^P^ 


6 


cluster  shovs  a set  of  practically  continuous  energy  levels  near  the  bandgap 
thus  obscuring  the  band-edge  structure  completely.  The  energy  levels  are 
much  more  numerous  than  what  the  reduction  of  the  S3nmnecry  suggests.  The 
fact  that  the  cluster  Is  less  symmetric  Is  not  a serious  problem  by  Itself 
since  we  can  relate  the  electron  states  of  the  cluster  to  those  of  the 
crystal  by  the  use  of  compatibility  relations  of  the  representations  between 
the  two  point  groups.  But,  for  any  cluster  model,  this  is  possible  only  for 
the  highest  symmetric  states. 

Now,  we  replace  all  of  the  six  outer  atoms  by  hydrogen  atoms  so  that 
they  do  not  have  any  dangling  bonds  of  their  own,  while  properly  containing 
the  bonds  originated  from  the  inner  phosphorus  and  gallium  atoms.  The  six 
electrons  of  the  hydrogen  atoms  are  expected  to  behave  similarly  as  the 
electrons  of  the  replaced  atoms  which  fill  the  bonds  directed  toward  the 
inner  phosphorus  or  gallium  atoms.  This  GaPH^  cluster  model  is  shown  in 
Fig.  1. 

B.  SCF  Xa  Multiple-Scattering  Calculation  for  the  GaPH^  Cluster 

The  sphere  radii  of  the  phosphorus  and  gallium  atoms  in  their 

O 

tetrahedral  configurations  are  not  much  different  from  each  ocher  (1.26  A 
for  Ga  and  1.10  A for  P).  Therefore,  equal  radii  were  used  for  both  the 
atomic  centers  of  the  GaPHg  cluster.  In  fact,  the  spherically  averaged 
potentials  of  the  ? and  Ga  atomic  centers  calculated  self-consistently  are 
equal  at  the  midpoint  between  their  nuclei.  The  radii  of  all  the  hydrogen 
atomic  centers  were  set  to  be  equal  to  the  same  radius  of  the  P or  Ga  atomic 
center  in  order  to  make  the  hydrogen  atoms  closely  related  to  the  atoms 
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^ they  replaced.  By  making  the  sizes  of  all  the  atomic  centers  equal, the 

Interatomic  region,  where  potential  Is  least  accurate  being  a constant.  Is 
kept  at  a minimum  for  a given  overlapping  scale  factor.  The  overlapping 
scale  factor  Is  the  ratio  of  the  overlapping  atomic  sphere  radii  to  the 
couching  atomic  sphere  radii  which,  in  this  case,  is  half  the  nearest- 
neighbor  distance.  This  factor  was  determined  by  satisfying  the  virlal 
theorem  (this  scheme  was  fully  explained  in  ref.  9).  For  the  GaPHg  cluster 
model  Che  overlapping  scale  factor  was  set  at  1.3  with  which  Che  virlal 

I 

theorem  was  satisfied  to  within  O.IZ. 

For  the  exchange  a parameters  Che  atomic  values  determined  by  Schwarz^ 

0.72620  and  0.70690,  were  used  for  the  P and  Ga  atomic  centers  respectively. 

18 

For  the  hydrogen  atomic  centers,  the  a value  of  0.77725  as  calculated  by 
the  spin-polarized  Xii  method  was  used  instead  of  the  Schwarz's  0.97304  which 
was  obtained  by  the  spin-restricted  Sx  method.  The  exchange  parameter  a 
for  the  interatomic  and  excramolecular  regions  of  Che  cluster  was  obtained 
by  averaging  the  a's  of  all  the  valence  electrons  from  the  phosphorus  and 
the  gallium  atoms;  it  is  0.71896  for  any  neutral  cluster  model  for  the  Ga? 
crystal.  The  electrons  of  the  hydrogen  atoms  were  not  accounted  for  since 
they  are  supposed  to  simulate  the  electrons  of  the  phosphorus  or  gallium 
atoms  in  the  cluster  model. 

The  basis  functions  were  constructed  from  the  s and  p functions  of 
the  P,  Ga,  and  H atoms.  In  Che  excramolecular  region  of  the  cluster,  the 
basis  functions  were  extended  to  have  s,  p,  d,  and  f type  functions. 


1 
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C.  Results  and  Discussions  ' 

The  electron  states  of  the  GaPHg  cluster  as  calculated  by  the  SCF  X& 
multiple-scattering  method  are  described  in  Table  I.  The  one-electron  energy 
levels  of  the  cluster  are  shown  in  Fig.  2(b),  compared  with  the  one-electron 
energy  levels  of  the  constituent  free  atoms  in  Fig.  2(a).  The  hydrogen's 
electron  levels  makeup  the  Rydberg  series,  and  the  electron  levels  of  the 
phosphorus  and  gallium  atoms  were  determined  by  the  SCF  Xa  method  using 
the  Schwarz's  a parameters.  Figure  2(c)  shows  the  energy-band  structure  of 
the  GaP  crystal  near  its  bandgap  (after  Dean^^) . 

In  the  calculation  for  the  GaPHg  cluster  it  was  found  that,  if  we 
fill  the  lowest  possible  electron  levels  by  the  available  14  valence  electrons, 
all  the  electron  sta'tes  of  the  top  valence  band  which  have  significant  con- 
tributions from  the  P 3p  and  Ga  4p  atomic  basis  states  are  unoccupied.  This 
resulted  in  a lack  of  proper  bonding  behavior  around  the  P and  Ga  atomic 
centers  and  over-population  of  electrons  in  the  Ga  and  H atomic  centers. 

The  reason  is  that  in  the  lower  valence  band  region  there  are  some  surface 
electron  states  which  have  their  electronic  charges  distributed  primarily 
in  the  boundary  hydrogen  atomic  centers.  They  do  not  correspond  to  any 
electron  states  of  the  GaP  crystal.  Therefore,  the  hydrogenic  surface 
electron  states  were  carefully  discerned  from  the  rest  of  the  valence  electron 
states  so  that  they  could  be  made  unoccupied  by  aiy  electrons  in  the  subsequent 
iterations.  The  le  and  2e  states,  whose  energy  levels  are  indicated  by 
dashed  lines  in  Fig.  2(b),  are  such  hydrogenic  surface  electron  states. 

After  they  were  made  unoccupied,  they  retained  their  characteristics  self- 
conslstently  in  the  rest  of  the  iterations. 
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Figures  3a  through  3k  show  the  contour  maps  of  the  wave  functions  of 
the  electron  states.  With  the  contour  maps  and  Table  I,  we  now  illustrate 
Che  electron  states  below  in  some  detail  in  the  order  of  increasing  eigen- 
energy  (refer  Fig  2(b)). 

The  la^^  state  is  mainly  made  of  a P 3s  basis  state  which  exhibits 
O bonds  with  the  Ga  4s  and  H-P  (H  atomic  centers  surrounding  the  P atomic 
center)  Is  basis  states.  The  Za^^  state  represents  a a bond  between  the 
H-Ga  Is  and  Ga  4s  basis  states.  There  is  also  a a bond  between  the  P 3s 
and  H-P  Is  basis  states, but  its  contribution  to  the  Za^^  state  is  relatively 
small. 

The  hydrogenic  surface  electron  state,  le,  is  almost  entirely  a 
H-Ga  Is  basis  state  containing  a weak  O bond  with  the  Ga  4p  basis  state. 

There  is  also  a very  weak  <7  bond  between  the  Ga  4p  and  P 3p  basis  states. 

The  other  hydrogenic  surface  electron  state,  2e,  is  mostly  a H-P  Is  basis 
state  showing  weak  a bonds  between  the  H-P  Is  and  P 3p  basis  states  and 
between  the  P 3p  and  Ga  4p  basis  states. 

For  the  la^^  state,  approximately  two  thirds  of  the  electronic  charge 
is  attributed  to  the  H-P  Is  basis  state  which  has  antibonding  behavior 
with  Che  nearby  basis  states.  It  carries  a moderate  amount  of  a 

a bond  between  the  Ga  4s  and  P 3s  basis  states.  The  3aj^  state  looks  like 
another  hydrogenic  surface  electron  state.  However,  it  has  less  contribution 
from  the  hydrogenic  basis  states  than  the  previously  determined  hydrogenic 
surface  electron  states,  le  and  2e.  Furthermore,  it  was  necessary  to  fill 
the  3a^  level  in  order  to  reproduce  the  features  of  the  band-edge  structure 
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of  Che  GaP  crystal  with  any  physically  reasonable  charge  distribution 
among  the  atomic  centers.  Therefore,  the  3a^  state  was  considered  to  be 
one  of  the  regular  valence  electron  states  of  the  cluster. 

Next,  Che  4a^  state  has  large  components  of  the  Ca  4s  and  H-Ga  Is 
basis  states  which  exhibit  an  antibond  between  them.  On  the  ocher  hand, 

Che  3e  state  represents  a strong  a bond  between  the  F 3p  and  Ga  4p  basis 
states,  and  the  5a^  state  represents  a strong  ir-bond  between  the  same  basis 
states.  This  set  of  the  three  highest  valence  electron  levels,  3e,  and 

of  Che  GaPHg  cluster  corresponds  to  the  top  valence  band  of  the  Ga? 
crystal.  They  all  have  a considerable  amount  of  their  electronic  charges 
in  the  interatomic  region  of  the  cluster.  Particularly,  the  top  valence 
electron  state,  has  more  than  one  third  of  its  electronic  charge  in 

the  interatomic  region. 

At  the  top  of  the  valence  band  of  the  GaP  crystal  there  is  a 
(or  t2)  state,  expressed  in  terms  of  a T^  single  group  representation,  which 
decomposes  into  an  a^^  and  an  e states  in  the  subgroup  symmetry.  There- 
fore, Che  Sa^^  state  of  the  cluster  has  the  correct  symmetry  as  the  top  of  the 
valence  band.  It  also  shows  the  correct  characteristics  of  a valence  band 


electron  state,  having  the  largest  contribution  from  the  P 3p 


basis 


state  with  a little  of  the  Ga  4p  and  forming  a strong  bond  between  them. 

Now  across  the  bandgap  there  are  several  excited  states.  The  first 
excited  state,  4e,  has  a relatively  large  component  of  the  extramolecular 
p basis  state,  It  puts  about  two  thirds  of  its  electronic  charge  in  the 


interatomic  region,  and  has  more  or  less  a uniform  electronic  charge  dis- 


u 


crlbuclon  over  Che  entire  space  of  the  cluster.  The  4e  state  shows  a 
weak  a bond  between  the  Ga  4p  and  H-Ga  2p' basis  states. 

The  6aj^  state  shows  sinilar  electronic  change  distribution  as  the 
4e  state  while  exhibiting  a a bond  between  the  Ga  4p  and  P 3p  basis  states. 
The  6aj^  state  has  the  sane  syranetry  as  the  state,  the  electron  state  of 
Che  crystal  at  the  boccon  of  Che  conduction  band  at  Che  center  of  Che 
Brillouin  zone.  Like  the  4e  state,  it  has  the  largest  contribution  from 
Che  Ga  4p  basis  state.  This  characteristic  is  the  sane  as  shown 

by  the  conduction  band  electron  states  of  Che  Ga?  crystal.  .Although  the  Sa^^ 
state  has  a higher  eigen-energy  than  the  4e  state,  it  is  aore  localized  in 
the  inner  ? and  Ga  atonic  centers.  The  ba^^  state  of  the  GaPH^  cluster 
represents  the  state  of  the  GaP  crystal,  and  the  4e  state  of  the  cluster 
corresponds  to  a state  near  t.he  state,  Che  electron  state  of  the  crystal 
at  the  absolute  bocton  of  the  conduction  band. 

Finally,  the  7aj^  state,  having  aore  chan  two  thirds  of  its  electronic 
charge  in  the  extranolecular  region,  shows  a very  strong  extraaolecular  s 
character.  This  state  has  also  sone  contributions  from  the  H-P  2s,  and 
H-Ga  2s  basis  states,  and  shows  a very  weak  bonding  for  all  the  aeons. 

The  7aj^  state  is  a true  surface  electron  state  of  the 

cluster. 

The  excited  states  have  their  wave  functions  more  concentrated  in 
the  interatoaic  and  extranolecular  regions  chan  the  valence  band  states, 
corresponding  to  the  sinilar  characteristics  in  the  crystal.  For  any 
cluster  model,  in  general,  the  higher  the  energy  level  of  a state  is,  the 


more  electronic  charge  it  has  in  the  interatomic  and  extramolecular  regions. 

Now  we  align  the  top  valence  electron  level,  5a. , of  the  cluster  to 
the  top  of  the  valence  band  of  the  crystal  and  linearly  expand  all  the 
other  levels  with  respect  to  this  by  the  factor  of  1.4  so  that  the  lowest 
excited  level,  4e,  of  the  cluster  is  also  aligned  to  the  bottom  of  the  con- 
duction band.  The  contraction  of  the  electron  energy-level  structure  of  the 

cluster  in  comparison  to  that  of  the  crystal  is  due  to  the  fact  that  there 
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are  only  limited  number  of  atoms  in  the  cluster.  The  linearly  expanded 
energy-level  structure  of  the  cluster  is  shown  in  rig.  2(d).  We  immediately 
see  in  the  Figs.  2(c)  and  2(d)  that  the  width  of  the  top  valence  band  of 
the  crystal  is  now  matched  very  closely  by  the  corresponding  band-width  of 
the  cluster  as  a consequence.  .A.C  the  same  time  the  oa^^  level  of  the  cluster 
also  fairly  well  matches  the  level  of  the  crystal.  These  correspond- 
ences of  the  levels  between  the  cluster  and  the  crystal  are  the  ones  that 
we  asserted  earlier  when  we  described  each  electron  state  of  the  cluster. 

Because  the  hydrogen  atoms  in  the  Ga?H,  cluster  represent  either 

Q 

phophorus  or  gallium  atoms  of  the  Ga?  crystal,  one  could  use  the  exchange 
parameter  a che  phosphorus  or  gallium  atoms  for  the  hydrogen  atomic 
centers.  This  change  of  the  a parameters  shifts  upward,  in  the  average  by 
0.025  Ry,  the  energy  levels  of  only  the  deep  valence  band  states,  through 
4aj^,  which  have  significant  components  of  the  hydrogenic  basis  states.  Since 
the  electron  states  near  the  bandgap  do  not  draw  any  significant  contributions 
from  such  hydrogenic  basis  states,  they  are  hardly  affected  by  the  change 


of  the  o values  in  the  hydrogen  atomic  centers. 


Consequently  the  interpretation 


of  ch«  physical  phenomena  associate  with  the  band-edge  structure 

will  be  the  same. 


1 

I The  total  energy  of  the  GaPH,  cluster  is  about  1.74  Ry  higher  than 

I O 

Che  sum  of  the  total  energies  of  the  constituent  free  atoms.  This  Indicates 
Chat  the  cluster  as  a molecule  is  unstable.  In  fact,  S*  uniform  expansion 

I 

of  the  GaPH,  cluster  decreases  the  total  energy  by  about  0.2  Ry  (in  magnitude) 

0 

I 

and  It  tends  to  be  monotonically  decreased  as  it  Is  further  expanded.  The 

unstableness  of  the  GaPH^  cluster  is  not  due  to  the  empty  hydrogenlc  surface 

0 

electron  states;  the  total  energy  of  the  cluster  with  filled  hydrogenlc 
surface  electron  states  Is  not  much  different  from  that  with  empty  hydrogenlc 
surface  electron  states. 

The  lonlcities  of  the  atomic  centers  shown  in  Table  II  were  obtained 
by  redistributing, for  each  valence  electron  state,  the  electronic  charge  in 
the  interatomic  and  extramolecular  regions  to  the  atomic  centers  in  pro- 
portion to  the  valence  electronic  charges  which  are  already  in  the  atomic 
centers  as  a result  of  SCF  calculation.  Note  that  all  the  hydrogen  atomic 
centers  have  practically  the  same  lonlcity.  If  we  make  the  hydrogens 
neutral  by  drawing  extra  electronic  charges  from  the  neighboring  P or  Ca  atomic 
center,  the  ionic ity  of  GaP  becomes  0.519  in  units  of  electronic  charge, 
j making  the  gallium  a cation  and  the  phosphorus  an  anion.  Since  all  three 

valence  electrons  of  the  gallium  atom  will  be  transferred  to  the  phosphorus 
atom  in  a completely  Ionic  model  of  GaP,  the  ionic  character  of  the  GaP 
crystal  may  be  said  to  be  about  17T  in  the  present  cluster  model  calculation. 

> 
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III.  NITROGEN  ISOELECTRON IC  ELECTRON  TRAP  IN  THE  CaP  CRYST.U, 

An  Isoeleccronlc  trap  is  an  isoelectronlc  aubstlcutlonal  atom  which 
induces  a bound  state  for  either  an  electron  or  a hole  in  the  bandgap  of 
a host  crystal.  Such  an  electron  trap  was  first  discovered  by  Thooas 
et  al.^^  for  nitrogen  substituents  at  phosphorus  sites  in  the  GaP  crystal. 
Nitrogen  iapurities  in  the  GaP  crystal  are  one  of  the  siaplest  kinds  of 
isoelectronic  traps.  Yet,  they  are  very  inportant  technologically  for  the 
efficient  luminescence  as  evidenced  by  the  green  llght-eniitcing  diode  made 
of  such  nitrogen-doped  GaP  crystal. 

The  effect  of  the  nicorgen  electron  crap  in  the  Ga?  crystal  can  be 
explained  qualitatively  as  follows.  A nitrogen  atom  attracts  an  electron 
more  strongly  chan  Che  phosphorus  atom  because  its  nucleus  is  less  screened. 
This  Is  reflected  in  the  high  electronegativity  of  the  nitrogen  atom  which 
is  3.0  compared  to  2.1  for  the  phosphorus  atom  in  Pauling's  scale.  Thus, 
the  nitrogen  impurity  produces  a bound  electron  state  which  is  highly 
localized.  Once  an  electron  is  crapped  in  the  state,  this  in  turn  gives 
rise  CO  a Coulomb  attractive  field  for  a weakly  bound  hole  whose  wave  function 
is  hydrogenic.  This  electron-hole  pair  conscitutes  a bound  excition.  Since 
Che  bound  electron  state  is  highly  localized  in  real  space,  its  k-vecCor 
components  are  widely  spread  in  Che  Brillouin  zone.  On  Che  ocher  hand,  Che 
hydrogenic  hole  state  is  very  much  localized  near  the  center  of  the  Brillouin 
zone  where  the  top  of  the  valence  band  of  the  crystal  is  located.  There- 
fore, once  Che  bound  exciCon  is  formed  due  to  the  nitrogen  impurity,  the 
electron  and  hole  can  directly  recombined,  emitting  a strong  photon  radiation 


15 


In  Che  green  part  of  Che  speccrum.  This  was  shown  In  an  opclcal  absorpclon 

and  fluorescence  experlaenc  as  A line  which  Is  abouc  0.0015  Ry  (0.02  eV) 

19 

below  Che  bandgap  energy. 

GaN  is  a dlrecc-gap  semlconduccor  which  cryscallizes  In  Che  wurczice 

sCruccure.  As  In  Che  zinc-blende  scruccure  of  GaP,  each  acom  in  GaN  is 

surrounded  by  four  CeCrahedrally  locaced  nearesc  neighbors  with  che  dis- 
■*  12 

Cance  of  1.94  A,  abouC  18Z  smaller  chan  chac  of  GaP.  Ic  is  inCerescing 

CO  noce  here  chac  che  bandgap  of  GaN  is  considerably  larger  (0.255  Ry  ac 

0®X)  Chan  chac  of  Ga?  (0.166  Ry  aC  0*I<)  , and  yec  a subscituced  nitrogen 

atom  causes  a bound  electron  state  in  the  bandgap  of  GaP. 

A heavier  isoelectronic  atom  of  smaller  electronegativity  chan  chat  of 

the  phosphorus  atom  could  produce  a localized'  bound  state  for  a hole  which 

Chen  attracts  an  hydrogenic  electron  by  Coulomb  interaction.  In  fact,  this 

20  21 

kind  of  exciton  was  found  in  bismuth-doped  GaP.  ’ In  this  case,  however, 
the  isoelectronic  impurity  does  not  lead  to  any  efficient  radiative  recom- 
bination since  the  localized  state  of  the  bound  hole  does  not  have  any 
significant  k-vector  components  near  the  X point  of  the  Brillouin  zone 
where  the  hydrogenic  electron  state  has  primary  k-vector  components.  Between 
phosphorus  and  bismuth  in  the  same  isoelectronic  column  of  che  periodic 
Cable,  there  are  arsenic  and  antimony  neither  of  which  has  been  found  to 

produce  any  bound  states  for  electrons  or  holes  when  substituted  at 

21 

phosphorus  site  in  the  GaP  crystal.  However,  their  electronegativities 
are  almost  equal  to  chat  of  che  bismuth  atom.  Therefore,  che  size  of  che 
substitutional  atom  must  be  another  important  factor  for  producing  a bound  state 

) 
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in  the  crystal. 

If  nitrogen  atoms  are  heavily  doped  in  GaP,  one  nitrogen  tends  to 

interact  with  another.  Depending  on  the  spacing  between  the  nitrogens, 

the  pairs  form  sharp  optical  peaks,  called  NN2>  etc.,  while  the  A 

line  caused  by  isolated  nitrogens  is  subsequently  weakened.  The  deepest 

line,  NNj^,  of  the  series  was  seen  at  0.0112  Ry  (0.152  eV)  below  the  bandgap 
19 

energy . 


A.  Calculational  Considerations 

In  order  to  study  the  nitrogen  electron  trap  in  GaP  we  replace 

the  phosphorus  atom  of  the  GaPH,  cluster  by  a nitrogen  atom.  This  makes  the 

0 

cluster  model  for  a nitrogen-doped  GaP  crystal  a GaNH^  cluster.  In 

• 

tetrehedtral  covalent  bonds  the  radius  of  the  nitrogen  atom  (0.70  A)  is 

* 16 

considerably  smaller  than  the  radius  of  the  phosphorus  atom  (1.10  A). 

However,  we  maintain  the  same  radii  for  all  the  atomic  centers  in  the 

GaNH,  cluster  so  as  to  be  consistent  with  the  GaPH-  cluster  model  for  the 
o 0 

GaP  crystal.  The  potential  of  the  N atomic  center  is  deeper  than  the 
potential  of  the  Ga  atomic  center  at  the  midpoint  between  them.  Therefore, 
if  we  had  used  the  scheme  of  determining  the  sizes  of  the  atomic  centers 
by  first  matching  the  potentials  at  the  point  of  contact,  the  size  of  the 
N atomic  center  would  have  been  even  larger  than  the  size  of  the  Ga  atomic 
center. 

The  exchange  parameter  a used  for  the  N atomic  center  is  0.75197  as 
calculated  by  Schwarz. But  the  average  a for  the  interatomic  and  extra- 
molecular  regions  of  the  GaPH^  cluster  is  kept  for  the  corresponding  regions 
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of  the  GaNH,  cluster  since  the  GaMH,  cluster  Is  still  supposed  to  simulate 
b b 

the  GaP  crystal  even  though  it  now  contains  an  impurity. 

It  is  a little  more  ambiguous  to  distinguish  the  hydrogenic  surface 

electron  states  from  the  regular  velance  electron  states  for  the  Gaini^ 

cluster  than  for  the  GaPH^  cluster.  However,  two  e states  of  the  GaHHg 

cluster  are  largely  the  same  in  their  xrave  function  characteristics  as 

the  two  e hydrogenic  surface  electron  states  of  the  GaPH^  cluster,  though 

their  eigen-energies  are  quite  different.  Since  the  wave  functions  of  the 

hydrogenic  surface  electron  states  would  not  be  modified  significantly  by 

the  presence  of  an  impurity  near  the  center  of  the  cluster,  we  identify 

the  two  e states  as  the  hydrogenic  surface  electron  states  of  the  GaNH^ 

cluster.  Such  close  relationships  maintained  between  the  GaNH  cluster  and 

0 

the  GaPHg  cluster  in  the  electronic  structure  as  well  as  in  the  geometric 
configuration  and  the  choice  of  exchange  parameters  ct's  will  differentiate 
the  GaNHg  cluster  from  a possible  cluster  model  for  a GaN  crystal. 

B.  Results  and  Discussions 

The  one-electron  energy  levels  of  the  GaNH^  cluster  are  shown  in 

Fig.  .4(b),  compared  with  those  of  the  constituent  free  atoms  of  the  cluster 

In  Fig.  4(a).  The  one-electron  energy  levels  of  the  GaPH^  cluster  are 

shown  again  in  Fig.  4(c)  for  convenience.  The  electron  states  of  the  GaHH, 

0 

cluster  are  described  in  Table  III,  and  the  contour  maps  of  their  ^rave 
functions  are  shown  in  Figs.  5a  through  5k.  In  the  following  we  will  study 
the  electron  states,  in  the  order  of  Increasing  eigen-energy,  while  establishing 
correspondence  relationships  between  them  and  the  electron  states  of  the 
GaPHg  cluster. 
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The  la^  state  Is  almost  entirely  an  N 2s  basis  state  which 

exhibits  weak  O bonds  with  the  Ga  4s,  H-N  Is,  and  H-Ga  Is  basis  states. 

This  state  corresponds  to  the  la^  state  of  the  GaPH^  cluster. 

The  2a state  has  most  of  its  electronic  charge  in  the  H-Ga  and 
Ga  atomic  centers,  thereby  making  a strong  a bond  between  the  H-Ga  Is 
and  Ga  4s  basis  states.  There  is  also  a weak  a bond  between  the  N 2s 
and  H-N  Is  basis  states.  This  2a.^  state  corresponds  to  the  2a^  state  of 
the  GaPHg  cluster. 

The  le  state  is  almost  entirely  a H-Ga  la  basis  state  containing 
weak  a bonds  between  the  H-Ga  Is  and  Ga  4p  basis  states  and  between  the 
Ga  4p  and  N 2p  basis  states.  Corresponding  to  the  le  state  of  the  GaPHg 
cluster,  it  is  a hydrogenic  surface  electron  state  of  the  GaNH^  cluster. 
Therefore,  this  state  was  made  unoccupied. 

The  3aj^  state  is,  in  a large  part,  an  antibonding  combination  of  the 
Ga  4s  and  H-Ga  Is  basis  states.  It  corresponds  to  the  4aj^  state  of  the 
GaPHg  cluster. 

The  2e  state  has  a strong  a bond  between  the  N 2p  and  Ga  4p  basis 

states,  and  a moderate  amount  of  a bond  between  the  N 2p  and  H-N  Is  basis 

states.  It  also  exhibits  a weak  T-bond  between  the  Ga  4p  and  H-Ga  2p  basis 

states.  This  2e  state  corresponds  to  the  3e  state  of  the  GaPH^  cluster. 

However,  it  has  more  contribution  from  the  H-N  Is  basis  state  than  the  3e 

state  of  the  GaPH,  cluster  does  from  the  corresponding  H-P  Is  basis  state. 

6 

Consequently  the  contribution  from  the  Ga  4p  basis  state  is  reduced  in  the 
case  of  the  GaNHg  cluster. 

The  4aj^  state  is  mostly  a H-N  Is  basis  state.  Therefore,  it  could 
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have  been  regarded  as  a hydrogenlc  surface  electron  state.  It  Is  instead 

treated  as  a regular  valence  electron  state  because  it  corresponds  to  the 

3a.  state  of  the  Ga?H,  cluster.  However,  the  contribution  of  the  N atomic 

center  for  the  4a^  state  comes  from  a hybrid  of  s and  p functions  unlike 

the  mainly  phosphorus  p function  for  the  Sa^^  state  of  the  GaPH^  cluster. 

Next,  the  3e  state  is  mostly  a H-N  Is  basis  state  which  i-uaws  a bonding 

behavior  with  the  N 2p  basis  state.  This  state  corresponds  to  the  2e 

hydrogenic  surface  electron  state  of  the  GaPH^  cluster;  it  is  the  other  of 

the  two  hydrogenic  surface  electron  states  of  the  GaNH^  cluster.  Unlike 

the  2e  state  of  the  GaPH,  cluster,  however,  the  3e  state  of  the  GaNH,  cluster 

0 0 

shows  a rather  strong  bond  between  the  N 2p  and  Ga  4p  basis  states,  and  the 

portion  of  the  N 2p  basis  function  which  is  responsible  for  the  bond  with 

the  H-N  Is  basis  function  is  the  inner  lobe,  i.e.  the  lobe  of  the  p function 

directed  toward  the  Ga  atomic  center. 

The  top  valence  state  5a^  corresponds  to  the  top  valence  state,  which 

is  also  5a,,  of  the  GaPH,  cluster.  As  such,  it  has  a moderate-to-strong 
± 0 

ir  bond  between  the  N 2p  and  Ga  4p  basis  states.  The  contribution  of  the 
N 2p  basis  state  to  the  5a,  state  of  the  GaNH,  cluster  is  larger  than  the 
contribution  of  the  corresponding  P 3p  basis  state  to  the  5aj^  state  of  the 
GaPHg  cluster. 

Now  across  the  bandgap,  the  6a^  state  represents  a strong  O bond 

between  the  Ga  4p  and  N 2p  basis  states,  corresponding  to  the  6a^  state  of 

the  GaPH,  cluster.  The  former  is,  however,  more  localized  chan  the  latter 
o 

inside  the  outersphere.  It  can  be  argued  chat  most  of  the  electronic  charge 
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In  the  excranolecular  region  and  some  in  the  P atomic  center  of  the  GaPH, 

0 

cluster  are  transferred  to  the  Ga  atomic  center  when  the  phosphorus  atom  is 

replaced  by  a nitrogen.  This  6a^  state  is  very  important  in  explaining  the 

isoelectronic  electron  trap  due  to  nitrogen  in  GaP  as  illustrated  later. 

The  4e  excited  state  of  the  GaNHg  cluster  corresponds  to  the  lowest 

excited  state  4e  of  the  GaPH,  cluster.  It  is  a little  more  localized  in 

0 

the  GalJHg  cluster  than  in  the  GaPH^  cluster  similarly  to  the  ba^^  state. 

This  4e  state  of  the  GaNHg  cluster  does  not  show  any  significant  bonding 

behavior  except  for  a weak  t bond  between  the  Ga  4p  and  H-Ga  2p  basis  states. 

Lastly,  the  Ta^^  excited  state  is  mostly  an  extramolecular  p state 

with  some  s character.  This  state  of  the  Ga^^Hg  cluster  will  correspond 

to  a state  of  the  GaPH,  cluster  whose  energy  level  is  slightly  higher  than 

0 

the  extramolecular  s type  Ta^^  level  of  the  GaPHg  cluster.  The  Ta^^  states 
are  true  surface  states  for  both  clusters. 

The  correspondence  relationships  of  the  electron  states  between  the 
GaNHg  cluster  and  the  GaPHg  cluster  are  indicated  by  the  interconnections 
between  Figs.  4(b)  and  4(c).  In  the  figures  the  electron  levels  of  the 
empty  hydrogenlc  surface  states  are  shown  as  dashed  lines  for  both  clusters. 
When  a nitrogen  atom  replaces  the  phosphorus  atom  in  the  GaPH^  cluster,  the 
energy  levels  are  in  general  shifted  down  because  the  potential  of  the 
nlgtrogen  is  deeper  than  the  potential  of  the  phosphorus  in  the  outer  region 
of  the  atomic  center.  However,  two  energy  levels,  4a^  and  3e,  of  the  GaNHg 
cluster  are  higher  than  the  corresponding  levels,  Sa^^  and  2e  respectively, 
of  the  GaPH,  cluster.  These  two  electron  states  of  the  GaNH,  clusteis  are 
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the  ones  which  have  large  partial-wave  components  of  the  H-N  Is  basis  states. 
The  increase  of  the  eigen-energies  results  from  a larger  Coulomb  repulsive 
interaction  in  the  H-N  atomic  centers  of  the  GaNH^  cluster  than  in  the 
corresponding  H-P  atomic  centers  of  the  GaPH^  cluster.  In  fact,  the  cal- 
culated ionicities  of  the  atomic  centers,  shown  in  Table  II,  indicate  a 
large  increase  of  the  electronic  charge  in  the  atomic  center  when  the 
nitrogen  atom  replaces  the  phophorus  atom  of  the  GaPH^  cluster. 

The  most  interesting  effect  of  the  nitrogen  substitutional  impurity 
is  to  shift  down  the  6a^  level  below  the  4e  level,  hence  into  the  bandga? 
region.  In  terms  of  the  energy-bands  of  the  crystal  it  can  be  explained  by 
the  following.  bound  electron  level  (the  6a^  level  of  the  GaNH^  cluster) 

' is  formed  a little  below  the  bottom  of  the  conduction  band  (the  4a  level 
of  the  GaPHg  cluster).  The  bound  electron  state  originates  from  the 
state  (the  6a,  state  of  the  GaPH,  cluster),  i.e.  the  bottom  conduction  band 
state  at  the  center  of  the  Brillouin  zone.  Unlike  the  usual  hydrogenic 
electron  states  due  to  donor  impurities,  the  bound  electron  state  due  to 
the  nitrogen  impurity  is  not  occupied  in  the  ground  state,  therefore,  it 
will  serve  as  an  electron  trap. 

As  discussed  earlier,  the  bound  electron  state  (6a^  of  GaMH^)  is 
more  localized  chan  the  original  state  (6a^  of  GaPH^) ; it  has  about  half 
of  its  electronic  charge  in  the  N atomic  center.  We  also  see  from  Table 
II  that  the  Ga  atomic  center  loses  0.763  electron  mainly  to  H-N  atomic  centers 
due  to  Che  presence  of  the  nitrogen  impurity.  This  indicates  that  the  Ga 
atomic  center  will  be  the  primary  site  for  an  extra  electron  which  can  be 
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crapped  by  the  nitrogen  impurity.  In  fact,  the  6a,  state  of  the  GaMH, 

1 0 

cluster  has  Che  largest  normalized  partial-wave  coefficient,  0.67,  of  a Ga 
4p  basis  state  among  all  the  electron  states  of  both  clusters. 

Once  an  electron  is  trapped,  a hole  will  be  attracted  by  the  negatively 
charged  nitrogen  impurity  complex,  thereby  making  a bound  exciton  which 
is  responsible  for  the  A line  in  optical  absorption  and  fluorescence  in  a 
nitrogen-doped  GaP  crystal. 


c 
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IV.  SUMMARY 

We  first  siaulated  a GaP  crystal  by  a GaPH^  cluster  model.  The 
SCF  Jfct  multiple-scattering  calculation  showed  that  the  electron  energy- 
level  structure  of  the  cluster  reproduces  the  essential  features  of  the  energy- 
band  structure  of  the  crystal.  The  electron  energy-level  structure  of  the 

GaPH,  cluster  was  expanded  linearly  to  fit  the  band-edges  of  the  GaP  crystal 

o 

in  order  to  compensate  the  lack  of  atomic  interaction  due  to  the  limited 

number  of  atoms  in  the  cluster.  It  is  remarkable  that  we  could,  at  least  ■ < 

qualitatively,  describe  the  crystal  by  using  such  a small  cluster  model  j 

(in  fact,  it  is  the  smallest-possible  neutral  cluster  model) . The  possible 
dangling  bonds  of  the  cluster  were  eliminated  by  using  the  hydrogen  atoms 
instead  of  three  phosphorus  and  three  gallium  boundary  atoms.  There  exist 
some  hydrogenic  surface  electron  states  in  this  cluster,  and  it  is  very  im- 
portant to  make  them  unoccupied  since  they  do  not  represent  any  electron 

states  of  the  crystal.  When  one  wants  to  work  with  a larger  cluster  j 

model  to  improve  quantitative  results,  the  use  of  boundary  hydrogen  atoms 
while  discerning  the  hydrogenic  surface  electron  states  carefully  will  still 
be  essential  in  order  to  get  correct  band-structure  of  the  crystal.  In  fact, 
this  scheme  should  be  used  for  cluster  models  for  any  covalent  crystal. 

The  isoelectronic  electron  trap  due  to  a nitrogen  Impurity  in  the 
GaP  crystal  was  studied  by  replacing  the  phosphorus  atom  of  the  GaPH^  cluster 
by  a nitrogen  atom,  thus  producing  a GaNH^  cluster, and  comparing  the  results 
of  the  SCF  Xa  multiple-scattering  calculations  for  the  two  clusters.  In  the 
GaMH,  cluster  a bound  electron  state  is  formed  in  the  bandgap,  about  0.05  Ry 
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below  Che  bottom  of  the  conduction  band.  It  is  obviously  impossible  in  the 
cluster  model  calculation  to  obtain  the  crap  level  with  such  an  accuracy  as 
can  be  found  in  an  experiment  which  suggests  the  level  is  about  0.0015  Ry 
below  the  bottom  of  the  conduction  band.  Because  the  bound  electron  state 
is  unoccupied,  it  acts  as  an  electron  trap.  Once  an  electron  is  trapped,  the 
nitrogen  impurity  complex  will  attract  a hole  by  Coulomb  interaction.  This 
electron-hole  pair  is  the  bound  exciton  which  causes  the  efficient  lumines- 
cence of  green  light  in  a nitrogen-doped  GaP  crystal,  as  indicated  by  the 
A line  in  optical  absorption  and  fluorescence  experiments. 

It  was  also  shown  that  the  bound  electron  state  originates  from  the 
state  which  corresponds  to  the  state  of  the  crystal.  In  the  crystal  the 
bound  electron  state  has  larger  k-vector  components  near  X point  than  near 
the  r point  in  the  Brillouin  zone.  However,  the  F point  contribution  of 
the  bound  electron  state  is  the  one  responsible  for  the  radiative  recombination 
of  the  bound  exciton.  In  this  sense,  the  cluster  model  successfully  explained 
the  luminescence  phenomena  of  the  nitrogen  impurity  in  the  GaP  crystal. 

The  GaNH,  cluster  model  is  different  from  a possible  cluster  model 
0 

for  a GaN  crystal  in  its  geometric  conf igviration,  particularly  in  the  lattice 
constant,  and  in  the  choice  of  the  exchange  parameter  a’s.  Also,  the  GaNH^ 
cluster  was  closely  related  to  the  GaPH^  cluster  in  determining  the  hydro- 
genic  surface  electron  states.  However,  it  seems  that  the  direct  bandgap 
nature  of  the  GaX  crystal  is  responsible  for  the  formation  of  the  bound 
electron  state  due  to  the  nitrogen  impurity  in  the  GaP  crystal. 
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TABLE  I 

Electron  states  of  the  CaPH^  cluster  as  calculated  by  the  SCF  Xa  multiple- 
scattering  method. 


They  are  listed  in  the  order  of  increasing  eigen-energy.  The  occupation 


numbers  of 

the  levels 

are 

shown  in  the  parentheses. 

Only  the 

Important  com- 

ponents  of 

the  normalized 

partial-wave  coefficients  and  elect: 

ronic  charge  dts- 

tributlons 

are  shown. 

(The  "In"  and  "Out"  represent 

the  interatomic 

and  extra- 

molecular 

regions  of 

the  1 

cluster  respectively.)  The 

constant 

potential  in 

the  interatomic  region  is 

-0.219  Ry. 

Electron 

Eigen- 

Partial-wave 

Elect 

:ronic 

charge 

state 

energy 

(Ry) 

coefficients 

distr 

■ibution 

la, (2) 

-0.379 

P 3s  H-P  Is  Ga  4s 

P 

H-P 

Ga 

0.81  0.39  -0.36 

0.65 

0.17 

0.13 

2a,  (2) 

-0.703 

H-Ca  Is  Ga  4s 

H-Ga 

Ga 

0.79  -0.42 

0.63 

0.13 

le  (0) 

-0.622 

H-Ga  Is 

H-Ga 

0.96 

0.92 

2e  (0) 

-0.572 

H-P  Is 

H-P 

( 

0.95 

0.90 

3a, (2) 

-0.511 

H-P  Is  P 3s 

H-P 

P 

1 

0.81  -0.33 

0.67 

0.12 

4a, (2) 

-0.351 

Ga  4s  H-Ga  Is  P 3s 

Ga 

H-Ga  P 

In 

L 

0.65  0.45  0.30 

0.43 

0.24  0 

.15  0.12 

3e  (4) 

-0.303 

P 3p  Ga  4p 

P 

Ga 

In 

0.75  -0.24 

0.57 

0.24 

0.15 

5a, (2) 

-0.226 

P 3p  Ga  4p 

P 

In 

Ga 

0.68  0.32 

0.46 

0.37 

0.09 

4e  (0) 

-0.103 

Out  p Ga  4p 

In 

Out 

Ga 

0.38  0.32 

0.67 

0.15 

0.10 

6a, (0) 

-0.066 

Ga  4p  P 3p  Out  p 

In 

Ga 

P Out 

0.51  0.41  0.32 

0.43 

0.26 

0.17  0.11 

7a, (0) 

-0.015 

Out  s H-P  2s  H-Ga  2s 

Out 

In 

H-P  H-Ga 

0.83  -0.16  -0.12 

0.71 

0.25 

0.024  0.015 
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TABLE  II 


loaicity  of  Che  atomic  centers  of  the  GaPH,  and  GaNH,  clusters  in 

o 0 

units  of  electronic  charge. 


GaPHg 

GaNHg 

Atomic  center 

lonicity 

Atomic  center 

lonicity 

P 

-1.170 

N 

-0.972 

Ga 

-0.136 

Ga 

0.627 

H-P 

0.217 

H-N 

-0.139 

H-Ga 

0.218 

H-Ga 

0.254 
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TABLE  III 

Electron  states  of  the  GaNH,  cluster  calculated  by  the  SCF  Xa  multiple- 

0 

scattering  method. 

The  are  listed  in  the  order  of  increasing  eigen-energy.  The  occupation 
numbers  of  the  levels  are  shown  in  the  parentheses.  Only  the  important  com- 
ponents of  the  normalized  partial-wave  coefficients  and  electronic  charge  dis- 
tributions are  shown.  (The  "In"  and  "Out"  represent  the  interatomic  and  extra- 
molecular  regions  of  the  cluster  respectively.)  The  constant  potential  in 
the  interatomic  region  is  -0.190  Ry. 


Electron 

state 

Eigen- 
energy  (Ry) 

Partial-wave 

coefficients 

Electronic  char 
distribution 

Se 

la^(2) 

-1.094 

N 2s 

0.96 

N 

0.92 

2a^(2) 

-0.844 

Xs  Gd  4s 

0.80  -0.54 

H-Ga 

0.65 

Ga 

0.29 

• 

le  (0) 

-0.736 

H-Ga  Is 

0.97 

H-Ga 

0.94 

3a  (2) 

-0.513 

Ga  4s  H-Ga  Is 

0.71  0.57 

Ga 

0.50 

H-Ga 

0.39 

2e  (4) 

-0.361 

N 2p  H-N  Is  Ga  4p 

. 0.70  -0.56  -0.34 

N 

0.49 

Ga 

0.12 

H-N 

0.11 

4aj^(2)  . 

-0.336 

H-N  Is  N 2p  N 3s 

0.87  0.20  0.18 

H-N 

0.77 

In 

0.15 

N 

0.07 

3e  (0) 

-0.315 

H-N  Is  N 2p  Ga  4p 

0.72  0.48  -0.32 

H-N 

0.53 

N 

0.23 

In 

0.14 

Ga 

0.10 

5a^(2) 

-0.276 

N 2p  Ga  4p  H-N  2s 

0.85  0.21  -0.19 

N 

0.72 

In 

0.17 

Ga 

0.05 

H-N 

0.05 

6a^(0) 

-0.142 

Ga  4p  N 2p 

0.67  -0.31 

Ga 

0.45 

In 

0.40 

N 

0.10 

4e  (0) 

-0.121 

Ga  4p  N 2p  Out  p H-Ga  2p 
0.43  0.32  0.27  0.21 

In 

0.57 

Ga 

0.19 

N 

0.10 

Out 

0.C8 

7a^(0) 

-0.011 

Out  p Out  s Ga  5p  H-Ga  2p 
0.71  -0.37  -0.23  0.17 

Out 

0.65 

In 

0.24 

Ga 

0.05 

H-Ga 

0.04 
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FIGURE  CAPTIONS 

Fig.  1.  The  cluster  model  of  GaPH^  for  a GaP  crystal.  The  small  spheres 
are  the  atomic  centers,  and  the  large  sphere  Is  the  outersphere 
outside  of  which  the  extramolecular  region  exists.  They  are  not 
shown  in  scale.  The  H-P  and  H-Ga  indicate  the  hydrogen  atomic 
centers  which  surround  the  P and  Ga  atomic  centers  respectively. 

Fig.  2.  One-electron  energy  levels  of  the  GaPH^  cluster  (a),  compared  with 
the  one-electron  energy  levels  of  the  constituent  free  atoms  of  the 
cluster  (b)  and  of  the  GaP  crystal  (after  Dean,  ref.  13)  (c).  The 
linearly  expanded  energy  levels  of  the  GaPH^  cluster  to  fit  the 
band-edges  of  the  GaP  crystal  are  shown  in  (d) . The  hydrogenic 
surface  electron  levels  of  the  cluster  are  shown  as  dashed  lines. 

represents  the  constant  potential  of  the  interatomic 

region  of  the  cluster. 

Fig.  3a.  Contour  map  of  the  la^^  wave  function  of  the  GaPKg  cluster.  It  is 
plotted  in  a plane  which  bisects  the  Ga  (middle  left) , the  P 
(middle  right),  a H-Ga  (upper  left),  and  a H-P  (lower  right)  atomic 
centers.  The  dashed  circles  are  the  projections  of  the  atomic 
centers  and  the  outersphere  of  the  cluster.  The  contour  curves 
were  selected  such  that  the  logarithm  of  the  magnitude  of  the  wave 
function  value  is  equally  spaced  for  the  consecutive  curves.  The 
maximum  wave  function  is  indicated  by  the  label  of  -fS  or  -5  de- 
pending on  the  sign  of  the  wave  function  value,  and  the  node  is 
indicated  by  the  label  of  either  +1  or  -1. 
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Fig.  3b.  Contour 

map 

of 

the 

2^1 

wave 

function 

of 

the 

GaPH, 

0 

cluster. 

Fig.  3c.  Contour 

map 

of 

the 

le 

wave 

function 

of 

the 

GaPHg 

cluster. 

Fig.  3d.  Contour 

map 

of 

the 

2e 

wave 

function 

of 

the 

GaPHg 

cluster. 

Fig.  3e.  Contour 

map 

of 

the 

2^1 

wave 

function 

of 

the 

GaPHg 

cluster. 

Fig.  3f.  Contour 

map 

of 

the 

4a^ 

wave 

function 

of 

the 

GaPHe 

cluster. 

Fig.  3g.  Contour 

map 

of 

the 

3e 

wave 

function 

of 

the 

GaPHg 

cluster. 

Fig.  3h.  Contour 

map 

of 

the 

wave 

function 

of 

the 

GaPH. 

0 

cluster. 

Fig.  3i.  Contour 

map 

of 

the 

4e 

wave 

function 

of 

the 

GaPHg 

cluster. 

Fig.  3j.  Contour 

nap 

of 

the 

6a^ 

wave 

function 

of 

the 

GaPH. 

0 

cluster. 

Fig.  3k,  Contour 

map 

of 

the 

wave 

function 

of 

the 

GaPH^ 

cluster . 

Fig.  4.  One-electron  energy  levels  of  the  GaNH,  cluster  (b),  compared  with 

o 

the  one-electron  energy  levels  of  the  constituent  free  atoms  of  the 
cluster  (a)  and  of  the  GaPH^  cluster  (c).  The  hydrogenic  surface 

electron  levels  are  indicated  by  dashed  lines.  represents 

the  constant  potential  in  the  interatomic  region  of  the  cluster. 

The  correspondence  relationships  of  the  electron  states  between  the 
GaNHg  cluster  and  the  GaPHg  cluster  are  indicated  by  the  inter- 
connections between  Figs,  (b)  and  (c). 

Fig.  5a.  Contour  map  of  the  la.  wave  function  of  the  GaNH,  cluster.  It  is 

L 0 

plotted  in  the  same  way  as  shown  in  Fig.  3a  except  that  the  N atomic 
center  now  replaces  the  P atomic  center. 


Fig.  5b.  Contour 

map 

of 

the 

2-1 

wave 

function 

of 

the 

GaNHg 

cluster. 

Fig.  5c.  Contour 

map- 

of 

the 

le 

wave 

function 

of 

the 

GaNHg 

cluster. 

Fig.  5d.  Contour 

map 

of 

the 

2-1 

wave 

function 

of 

the 

Ga.NHg 

cluster. 

Fig.  5e.  Contour 

map 

of 

the 

2e 

wave 

function 

of 

the 

GaNHg 

cluster. 
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Fig.  5f.  Contour 

map 

of 

the 

4ai 

wave 

function 

of 

the 

GaNHg 

cluster. 

Fig.  Sg.  Contour 

map 

of 

the 

3e 

wave 

function 

of 

the 

GaNH, 

D 

cluster. 

Fig.  5h.  Contour 

map 

of 

the 

wave 

function 

of 

the 

GaNHg 

cluster. 

Fig.  51.  Contour 

map 

of 

the 

6a^ 

wave 

function 

of 

the 

GaNHg 

cluster. 

Fig.  5j.  Contour 

map 

of 

the 

4e 

wave 

function 

of 

the 

GdNHg 

cluster. 

Fig.  5k.  Contour 

map 

of 

the 

wave 

function 

of 

the 

GaNHg 

cluster . 

